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METHOD AND SYSTEM FOR PREDICTING NUCLEIC 
ACID HYBRIDIZATION THERMODYNAMICS AND 
COMPUTER-READABLE STORAGE MEDIUM FOR USE THEREIN 



CROSS-REFERENCE TO RELATED APPLICATION 



5 This application claims the benefit of U.S. provisional application 

Serial No. 60/209,778, filed June 7, 2000, entitled "Nucleic Acid Hybridization 
Prediction and Biochemical Techniques Utilizing Same." 



STATEMENT REGARDING FEDERALLY SPONSORED 
RESEARCH OR DEVELOPMENT 



10 This invention was made with government support under Grant No. 

HG02020 provided by NIH. The United States Government has certain rights in 
the invention. 



BACKGROUND OF THE INVENTION 



1 . Field of the Invention 



15 This invention relates to methods and systems for predicting nucleic 

acid hybridization thermodynamics and computer-readable storage medium for use 
therein. 



2. Background Art 

Improvement of the efficiency of hybridization-based techniques 
20 requires the optimization of the binding between two sequences. Accurate 
prediction of the thermodynamics allows optimal choice of the sequences, 
temperature, and salt conditions. Hence, the prediction of nucleic acid 
thermodynamics is important to optimize techniques like PGR (Saiki et al., 1988), 
Southern and Northern blotting (Southern, 1975), antigene targeting (Freier, 1993), 
25 and Kunkel site-directed mutagenesis (Kunkel et aL, 1987). 
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Hybridization prediction is also important for designing DNA 
microchips that have a wide field of application ranging from diagnostics (Hacia, 
1999; Yershov et al., 1996) to gene expression analysis (Ferea et al., 1999) and 
drug discovery (Debouk and Goodfellow, 1999). Microchips contain a large 
5 number of DNA probe sequences that have to be designed to specifically hybridize 
target sequences in a pool of DNA fragments. First, a DNA probe should be 
designed to bind to only one site of only one DNA target. Second, the different 
DNA probe sequences need to hybridize to their targets under the same temperature 
and solution conditions. Moreover, in sequencing by hybridization (Fodor et al., 
10 1993; Mirzabekov, 1994) where microchips are used to determine the sequence of 
given DNA, one has to be able to know hybridization thermodynamics to 
discriminate signals resulting from perfectly matched and mismatched probe/target 
hybridizations. 

Another widely used technique that requires hybridization prediction 
15 is the fluorescence in situ hybridization (FISH) technique (Gall and Pardue, 1969). 
In this technique, a fluorescent tagged nucleic acid probe is designed to specifically 
hybridize cellular or tissue section nucleic acids. The target of these probes can 
either be endogenous DNA, messenger RNA or viral and bacterial sequences. 

Therefore, FISH is used to monitor gene expression (McNicol and 
20 Farquharson, 1997), detect infectious agents (Bashir et al., 1994; McNicol and 
Farquharson, 1997; Pollanen et al., 1993), study cell cycle (McNicol and 
Farquharson, 1997), map chromosomes and study nuclear architecture (Heng et al., 
1997). It was also determined that a set of probes can be used simultaneously 
(multiFISH) to detect different loci (Pagon, 1997). Once again, prediction of 
25 hybridization is essential to insure specificity. Nucleic-acid hybridization prediction 
is also important for the design of oligonucleotide aptamers or antisense 
oligonucleotides (Cohen, 1992) that can be used for various therapeutic applications. 
A new type of probes known as molecular beacons (Bonnet et al., 1999; Tyagi et 
al. 1998) that are very specific has been developed and shown to be efficient for 
30 mutation analysis (Giensendorf et al., 1998) and multiplex detection of single 
nucleotide variations (Marras et aL, 1999). The design and prediction of the 
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thermodynamics of these beacons is helped by hybridization thermodynamics 
prediction (Bonnet et al., 1999). Accurate prediction of hybridization is also 
important for the practical realization of DN A-based or more generally nucleic acid- 
based computers. (Adleman, L.M., 1994). 

5 The development of molecular biology techniques based on 

hybridization (PCR, FISH, DNA microchips, etc.) has resulted in a need for 
efficient automated ways to design probes and primers. In the last decade, 
numerous algorithms have been developed to optimize the design of primers and 
probes for various applications (Rychlik and Rhoads, 1989)(Breslauer et al., 1986; 

10 Chen and Zhu, 1997; Dopazo et al., 1993; Haas et al., 1998; Hillier and Green, 
1991; HyndmanetaL, 1996; LietaL, 1997; Linket al., 1997; Pesoleet aL, 1998; 
Proutski and Holmes, 1996). Numerous unpublished software to predict primers 
are also made available by research groups and biotech companies on the World 
Wide Web (Primer3 from the Whitehead Institute for Biomedical Research. Primer 

15 Express™ from PE Biosystems, DNAstar from IDT, etc.) 

There are currently many software packages on the market for DNA 
primer design including: OLIGO, PRIMER PREMIER, OSP, GCG, PrimerMaster, 
and Primo. None of the current programs, however, were written by experts in 
DNA thermodynamics; thus, there are many improvements that can be made. 

20 Nearly all of the current software packages contain mistakes that result from a lack 
of understanding of the underlying theory of DNA hybridization. PCR is a fairly 
robust process and thus even crude programs make predictions that work 90-95 % 
of the time. Multiplex PCR primer design, however, is not at all trivial and 
detailed knowledge of the physical chemistry of DNA hybridization, and the 

25 availability of an accurate thermodynamic database are essential to reliable design 
of multiplex PCR primers. In multiplex PCR, several primers must be designed to 
specifically bind to different sites on target DNA at a given temperature with 
minimal background binding to mismatch sites and with minimal cross- 
hybridizations between pairs of primers. The design of molecular beacons for DNA 

30 oligonucleotide arrays is also very challenging because of the complex competing 
equilibria. 
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Most of the existing programs aiming at finding an optimum probe 
that binds a specific location on a target, however, do not include accurate stability 
rules for hybridization and neglect or poorly approximate competitive binding sites, 
strand folding and strand dimerization. 

5 U.S. Patent Nos. 5,593,834 and 6,027,884 to Lane et al. disclose 

methods to design and construct DNA sequences with selected reaction attributes. 

In summary, prediction of nucleic acids thermodynamics is important 
to optimize various molecular biology techniques including multiplex PCR, DNA 
microchips, molecular beacons, and fluorescence in situ hybridization. Most of the 
10 available programs for probe design do not include a complete parameterization and 
often do not account for mismatches. Moreover, single strand folding is not taken 
into account, which often leads to inaccurate predictions. 

SUMMARY OF THE INVENTION 

An object of the invention is to provide a method and system for 
15 predicting nucleic acid hybridization thermodynamics and computer-readable storage 
medium for use therein wherein the invention utilizes a thermodynamically rigorous 
approach to evaluate the quality of probes and simulate probe/target hybridization. 

Another object of the invention is to provide a method and system for 
predicting nucleic acid hybridization thermodynamics and computer-readable storage 
20 medium for use therein wherein the invention also takes into account single strand 
folding thermodynamics to calculate effective hybridization thermodynamics. 

In carrying out the above objects and other objects of the present 
invention, a method for predicting nucleic acid hybridization thermodynamics is 
provided. The method includes providing a database of thermodynamic parameters, 
25 receiving hybridization information which represents at least one sequence, 
receiving correction data, receiving a first set of data which represents hybridization 
conditions, and calculating hybridization thermodynamics including net 
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hybridization thermodynamics based on the hybridization information, the 
thermodynamic parameters, the correction data and the first set of data. 

The hybridization thermodynamics of individual single stranded, 
bimolecular and higher order complexes may be statistically weighted in a numerical 
5 process and the equilibrium concentration of each species is output. 

The correction data may include folding correction data and/or linear 
correction data. 

The thermodynamic parameters may include DNA thermodynamic 

parameters. 

10 The DNA thermodynamic parameters may include dangling end 

parameters and/or coaxial stacking parameters. 

The DNA thermodynamic parameters may further include terminal 
mismatch parameters. 

The thermodynamic parameters may include RNA thermodynamic 
15 parameters and/or hybrid DNA/RNA thermodynamic parameters. 

The thermodynamic parameters may further include DNA loop 
thermodynamic parameters. 



20 calculated for the duplex. 

The hybridization information may further represent at least a section 
of a target and a length of at least one primer or probe complimentary to the target. 



The hybridization information may represent top and bottom strand 
sequences which form a duplex and wherein the hybridization thermodynamics are 
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The hybridization thermodynamics may be calculated for a plurality 
of primers or probes complimentary to the target. 

The hybridization information may represents at least a section of a 
target and a primer or probe. 

5 A length of the target may be longer than a length of the primer or 

probe and wherein the hybridization thermodynamics are calculated for a best 
target/primer or target/probe complex and for competitive mismatch complexes. 

Hybridization information may represent at least a section of a target 
and a primer or probe and wherein a length of a target is longer than the length of 
10 the primer or probe and wherein the hybridization thermodynamics are calculated 
for a best target/primer or target/probe complex and for competitive target/primer 
or target/probe complexes. 

The method may further include calculating concentration of each 
species in a solution at a plurality of temperatures. 

15 Hybridization information may also represent a primer or probe and 

wherein the length of the target is longer than a length of the primer or probe and 
^ wherein the hybridization thermodynamics are calculated for a best target/primer 
or target/probe complex and for competitive mismatch complexes and wherein the 
method may further comprise calculating concentration of every species in a 

20 solution at a plurality of temperatures. 

The hybridization thermodynamics may be calculated for at least two 
best target/primer or target/probe complexes and for their corresponding 
competitive mismatch complexes and wherein the method may further comprise 
correcting for any interactions between the at least two best target/primer or 
25 target/probe complexes and their components. 
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Further in carrying out the above objects and other objects of the 
present invention, a system for predicting nucleic acid hybridization 
thermodynamics is provided. The system includes a database of thermodynamics 
parameters, means for receiving hybridization information which represents at least 
5 one sequence, and means for receiving correction data. The system further includes 
receiving a first set of data which represents hybridization conditions, and means 
for calculating hybridization thermodynamics including net hybridization 
thermodynamics based on the hybridization information, the thermodynamic 
parameters, the correction data and the first set of data. 

10 The hybridization thermodynamics of individual single stranded, 

bimolecular and higher order complexes may be statistically weighted in a numerical 
process and the equilibrium concentration of each species is output. 

The correction data may include folding correction data and/or linear 
correction data. 



15 The thermodynamic parameters may include DNA thermodynamic 

parameters such as dangling end parameters. 

The DNA thermodynamic parameters may include coaxial stacking 
parameters and/or terminal mismatch parameters. 

The thermodynamic parameters may include RNA thermodynamic 
20 parameters and/or hybrid DNA/RNA thermodynamic parameters. 

The thermodynamic parameters may further include DNA loop 
thermodynamic parameters. 

The hybridization information may represent top and bottom strand 
sequences which form a duplex and wherein the hybridization thermodynamics are 
25 calculated for the duplex. 
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The hybridization information may also represent at least a section 
of a target and a length of at least one primer or probe complimentary to the target. 

The hybridization thermodynamics may be calculated for a plurality 
of primers or probes complimentary to the target. 

5 The hybridization information may represent at least a section of a 

target and a primer or probe. 

A length of the target may be longer than a length of the primer or 
probe and wherein the hybridization thermodynamics are calculated for a best 
target/primer or target/probe complex and for competitive mismatch complexes. 

10 Hybridization information may represent at least a section of a target 

and a primer or probe and wherein a length of a target is longer than the length of 
the primer or probe and wherein the hybridization thermodynamics are calculated 
for a best target/primer or target/probe complex and for competitive target/primer 
or target/probe complexes. 

15 The system may further include means for calculating concentration 

of each species in a solution at a plurality of temperatures. 

Hybridization information may also represent a primer or probe and 
wherein the length of the target is longer than a length of the primer or probe and 
wherein the hybridization thermodynamics are calculated for a best target/primer 
20 or target/probe complex and for competitive mismatch complexes and wherein the 
system may further comprise means for calculating concentration of every species 
in a solution at a plurality of temperatures. 

The hybridization thermodynamics may be calculated for at least two 
best target/primer or target/probe complexes and for their corresponding 
25 competitive mismatch complexes and wherein the system may further comprise 
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means for correcting for any interactions between the at least two best target/primer 
or target/probe complexes and their components. 

Still further in carrying out the above objects and other objects of the 
present invention, a computer-readable storage medium having stored therein a 

5 database of thermodynamics parameters and a computer program are provided. The 
computer program executes the steps of: a) receiving hybridization information 
which represents at least one sequence; b) receiving correction data; c) receiving a 
first set of data which represents hybridization conditions; and d) calculating 
hybridization thermodynamics based including net hybridization thermodynamics 

10 based on the hybridization information, the thermodynamic parameters, the 
correction data and the first set of data. 

The hybridization thermodynamics of individual single stranded, 
bimolecular and higher order complexes may be statistically weighted in a numerical 
process and the equilibrium concentration of each species is output. 

15 The correction data may include folding correction data and/or linear 

correction data. 

The thermodynamic parameters may include DNA thermodynamic 

parameters. 

The DNA thermodynamic parameters may include dangling end 
20 parameters and/or coaxial stacking parameters. 

The DNA thermodynamic parameters may further include terminal 
mismatch parameters. 

The thermodynamic parameters may include RNA thermodynamic 
parameters and/or hybrid DNA/RNA thermodynamic parameters. 
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The thermodynamic parameters may further include DNA loop 
thermodynamic parameters. 

The hybridization information may represent top and bottom strand 
sequences which form a duplex and wherein the hybridization thermodynamics are 
5 calculated for the duplex. 

The hybridization information may represent at least a section of a 
target and a length of at least one primer or probe complimentary to the target. 

The hybridization thermodynamics may be calculated for a plurality 
of primers or probes complimentary to the target. 

10 The hybridization information may represent at least a section of a 

target and a primer or probe, 

A length of the target may be longer than a length of the primer or 
probe and wherein the hybridization thermodynamics are calculated for a best 
target/primer or target/probe complex and for competitive mismatch complexes. 

15 Hybridization information may represent at least a section of a target 

and a primer or probe and wherein a length of a target is longer than the length of 
the primer or probe and wherein the hybridization thermodynamics are calculated 
for a best target/primer or target/probe complex and for competitive target/primer 
or target/probe complexes. 

20 The program may further execute the step of calculating 

concentration of each species in a solution at a plurality of temperatures. 

Hybridization information may also represent a primer or probe and 
wherein the length of the target is longer than a length of the primer or probe and 
wherein the hybridization thermodynamics are calculated for a best target/primer 
25 or target/probe complex and for competitive mismatch complexes and wherein the 
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program may execute the step of calculating concentration of every species in a 
solution at a plurality of temperatures. 

The hybridization thermodynamics may be calculated for at least two 
best target/primer or target/probe complexes and for their corresponding 
5 competitive mismatch complexes and wherein the program may execute the step of 
correcting for any interactions between the at least two best target/primer or 
target/probe complexes and their components. 

The above objects and other objects, features, and advantages of the 
present invention are readily apparent from the following detailed description of the 
10 best mode for carrying out the invention when taken in connection with the 
accompanying drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIGURE 1 is a schematic drawing wherein multiple equilibria are 
considered for concentration calculations; 

15 FIGURE 2a is a schematic drawing of a user input interface wherein 

the user provides various input information for a first module of the invention; 

FIGURE 2b is a schematic drawing of a user output interface wherein 
a computer provides output information corresponding to the input information of 
Figure 2a; 

20 FIGURE 3a is a schematic drawing of a user input interface wherein 

the user provides various input information for a second module of the invention; 

FIGURE 3b is a schematic drawing of a user output interface wherein 
a computer provides output information corresponding to the input information of 
Figure 3a; 
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FIGURE 4a is a schematic drawing of a user input interface wherein 
the user provides various input information for a third module of the invention; 

FIGURE 4b is a schematic drawing of a user output interface wherein 
a computer provides output information corresponding to the input information of 
5 Figure 4a; 

FIGURE 5a is a schematic drawing of a user input interface wherein 
the user provides various input information for a fifth module of the invention; 

FIGURE 5b is a schematic drawing of a user output interface wherein 
a computer provides output information corresponding to the input information of 
10 Figure 5a; 

FIGURE 6 is a block diagram flow chart illustrating the solution of 
conservation equations of the present invention; 



FIGURE 8 shows prediction of molecular beacon net hybridization 
1 5 thermodynamics ; 

FIGURE 9 shows simulation of molecular beacon hybridization 
concentrations at temperatures from 0 to 100°C; 



FIGURE 11 shows match vs. mismatch hybridization simulation at 
20 different temperatures; 

FIGURE 12 shows a general case of competitive hybridization 
equilibria that can be solved using the described numerical methods; and 



FIGURE 7 is a schematic diagram illustrating multiplex PCR design; 



FIGURE 10 is a diagram of match vs. mismatch hybridization; 
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FIGURE 13 is an example of simultaneous equations for the general 
five molecule case. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

In general, the method and system of the present invention include 
5 rigorous thermodynamic parameterization for Watson-Crick base pairs, internal 
mismatches, terminal mismatches, terminal dangling ends, co-axial stacking 
interactions, sodium and magnesium salt dependence, denaturants (urea, formamide, 
DMSO). In addition, loop parameters for hairpins, internal loops, bulges, and 
multibranched loops are included. For DNA essentially all the parameters have 

10 been previously published or all included in the Appendix hereto. Specifically, the 
parameters which have been published include Watson-Crick parameters, sodium 
dependence, GT, GA, CT, AC, AA, CC, GG, and TT mismatches. The 
parameters included herein include dangling ends, terminal mismatches, DNA loop 
parameters, and co-axial stacking parameters. For RNA, the parameters have been 

15 published by Douglas H. Turner et al. For DNA/RNA hybrid duplexes, the 
parameters have been published by Naoki Sugimoto. 

The method and system are adapted for future implementation of 
parameters for modified nucleosides (including but not limited to inosine, 5- 
nitroindole, PNA, MOE-modified RNA, and iso-bases). With these parameters, it 

20 is possible to predict the melting temperature, Tm, of a duplex within 2°C on 
average. Correction for surface effects for DNA chip arrays is also implemented. 
In addition to predicting duplex hybridization, the software accounts for single- 
strand secondary structure. This is accomplished by a new numerical procedure for 
solving complex coupled equilibria (multi-state model). With this approach, it is 

25 possible to accurately predict not only the Tm for hybridization but also the 
concentration of every species in the solution (e.g. match duplex, mismatch, duplex, 
folded target, folded primer, primer dimer, etc.) at every temperature from 0 to 
100°C. Thus, it is possible to use this software to design oligonucleotide 
hybridization with optimized temperature, salt, and strand concentrations. 
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Predicting Accurately Primer Target Interaction Stability 

The stability of a primer/target or probe/target complex can be 
described by the free energy of association of the probe and the target. The most 
accurate way to calculate free energy of association is to use the nearest-neighbor 
5 model with accurate thermodynamic parameters. Thermodynamic parameters 
should account for Watson-Crick base pairs (SantaLucia et al., 1996; Allawi & 
SantaLucia, 1997; SantaLucia, 1998), single mismatches (Allawi & SantaLucia, 
1997), terminal mismatches (disclosed herein), dangling ends (Bommarito, Pugret 
& SantaLucia, 2000) and possibly double mismatches. Proper calculation of the 

10 monovalent and divalent salt dependence is also important (SantaLucia, 1998). 
Other loop motifs for hairpins, bulge, internal loops and multi-branched loops are 
important for single strand secondary structure prediction, but are often very 
crudely approximated. Moreover, when primer and target folding can occur, a set 
of coupled equilibria should be used to model the system. The nearest-neighbor 

15 model needs to be used to determine the equilibrium constant of each equilibrium. 
The determination of possible primer or target folding can be addressed by using 
secondary-structure prediction algorithms like M. Zuker's MFOLD (Zuker, 1989). 

Secondary Structure and Net Hybridization Thermodynamics 

Species Concentration Calculations 

20 Consider a system of strands SI and S2 with four states: folded 

target, folded probe, probe bound to target, and random coil target and probe. The 
model can be described by three equilibria as shown in Figure 1. 

The concentrations of every species for such a system can be 
analytically determined. The three equilibrium constants for such a system are 
25 shown below: 

S1+S2-DH K, = [DH] (1) 

1 [SI] [52] v } 
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S1SH1 K 2 - ^ (2) 



S2=H2 * 3 = ^ (3) 

where SI, S2, HI, H2, and DH are the random coil SI, the random coil S2, the 
folded strand HI, the folded strand H2, and the double helix DH, respectively. The 
5 conservation of SI and S2 leads to the following equations: 

C T s °' al = Sl+Hl+DH (4) 

eg" 1 = S2+ H2+ DH (5) 

where Csf 1 are the total concentrations of SI and S2. [DH] and [S2] can be 
10 expressed as a function of [SI] by substituting the [HI] obtained from Equation 2, 
in Equation 6, and then substituting the [DH] obtained by Equation 6 in Equation 
1. 



[DH]=C™-[S\-]-K 2 [Sl] (6) 



[«] = ^ =, [S2] - 4r-rtt]-*CT (7) 



15 Substitution of [H2], [DH], and [S2] from Equations 3, 6 and 7 in 

Equation 5 leads to an expression of Ki that can be rearranged as a quadratic 
equation in [SI]: 

[SI] 2 (K, + K 2 K L ) + [S1](K, C S2 + K 3 + K 2 K 3 - K, C S1 + K 2 + 1) - (K 3 + 1) C S1 = 0 

(8) 

20 This equation is simplified by making the following substitutions: 
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a = (K, + K 2 K.) (9) 
b = (K t C^ 2 0tal + K 3 + K 2 K 3 - K, Cj° tal + K 2 + 1) (10) 
c = (K 3 + 1) C s T r' (ID 



The physical solution of the quadratic equation (i.e. positive root) is (Press, 1999): 



5 [51]= -^f^ (12) 

2a 



or 



lS1] -- t -Z-4<,c (13> 



The second equation has better numerical stability (Press, 1999). 
[DH], [S2], [HI], and [H2] can then be calculated using Equations 1-3. 
10 Determination of Net Free Energy 

The net free energy of hybridization is calculated as follows: 



= -RT\nK nel (14) 



where 



[DH] 

single stranded -» L single stranded J 



15 where [SI single stranded] and [S2 single stranded] are the concentrations of SI and 
S2 either in the random coil state or the hairpin states, at the temperature of the 
simulation. Using the conservation of SI and S2, Equation 14 is rewritten as 
follows: 
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&G% net = RT ^ {C r^_ [D ^ ) ^T^_ [Dm) (16) 

Note that AG° rm has the unusual property that it depends on the total strand 
concentrations, and Cg* . The net free energy expresses the duplex formation 
equilibrium free energy corrected for secondary-structure formation in the single 
strands. 



Determination of Net Melting Temperature 

If the strands are non self-complementary two cases have to be 
considered depending on the relative strand concentrations: 

1) If SI is the limiting reagent (Cfg* < C T s T l ), at T M : 



10 [DH]=~C^ aI (17) 

The concentrations of strands [SI] and [S2] are given by the following relations: 



Cir 1 = \clT l + K 2 [SI] + [51] (18a) 



r 



Total 



Cjr* = !C"'+A: s [S2]+[S2] (19a) 



(i Total _ }_{-iTotal 

LJ [52]= * (19b) 

K 3 + 1 
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The replacement of [SI] and [S2] in Equation 1 gives: 



[SI] [52] r-iTotal 1 Total 

2 



o - \r Total r Tot ° l _l Kl + K i K i + ^ 3 + 1 (20b) 



Using the relation AG° T = -R T In K, Equation 20 is arranged as follows: 



1 f ~AGy(2) -AG° r (2)-AG° r (3) -AG° r (3 ) \ AG<VO) 



0= ^-Cj°""- Cjf + 1 e «• + e «■ + e «• + i\ e ~ «* 



(21) 

AG° T can then be decomposed as AG° T =AH°- T AS° (assuming AC° 0 ) to obtain: 



1 -Afl°(2)+Aff°(l) +/VS°(2)-AS°(1) 

TVi*^/ T~ + „1 : 

2 



-M°(2)-A//°(3)fAff'(l) AS«(2)+Ag«(3)-AS»0) -A#°(3)+A//°(l) A5°(3)-A5°(l) A/T»(l) -AS»(1) 

(22) 

10 The above equation can be solved by bisection or other numerical techniques to find 
T. This solution is the net melting temperature. 

2) If S2 is the limiting reagent (C^f 1 < C^ tal ), the following relation can 
be deduced by a similar approach: 



n - 1 r Tota{ r Total j. K * + K * K * + K * + 1 

15 Again, application of the bisection method to an equation symmetric to Equation 22 
affords the net melting temperature. 
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If the strand S is self complementary, the reactions are described by 
the following equilibria: 



S + S*DH ^ = ^ (24) 



S*H * 2 = ^ (25) 



AtT M :[DH]=\cl°' a ' (26) 
The strand conservation equation is: 

C™ = [S]+[H]+2[DH] (27) 
Insertion of [H] and [DH] from Equations 25 and 26 in Equation 27 leads to: 



si Total 

Cf'"' = 2(K 2 + o [S] = 2( ^ + i) (28) 



10 Introduction of [S] in Equation 24 gives: 



Kl = W = Cf" ° 0 = + 2 ^2 " ^Cf + 1 (29) 



Using the relation AG° T = -R T In K, Equation 29 is rearranged as follows: 



-2AG° r (2) -AG' r (2) AG' r (l) 

0 = e « r + 2e * r - C T s °' al e RT + 1 (30) 
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AG° T can then be decomposed as AG° T =AH°- T AS° to obtain: 

-2A/f°(2) +2A5°(2) -&H°(2) +A5°(2) -AH°{\) +AS°(1) 

0=e RT e R +2e RT e R -C T s otal e RT e R +1 (31) 

This equation can be solved by bisection to afford the net melting temperature. 

An experimentally validated example of the accuracy of the net 
5 hybridization thermodynamics is shown in Figure 8 for molecular beacons. At the 
top of Figure 8 are the predicted thermodynamics for simple duplex formation 
assuming no competing single strand secondary structure. Using Module 1 of the 
invention, these results are similar to what would be predicted using other 
commercial software (such as oligo 6.0), though our thermodynamic database 

10 includes the dangling end effects and salt corrections are more accurate than other 
software. The middle of Figure 8 shows the single strand folding at the molecular 
beacon as output from DNA-MFOLD. The bottom table of Figure 8 shows the 
experimentally determined A6 (effective) and Tm (effective) published in Bonnet et 
al. 1999, as well as the effective Tm and A6 (effective) predicted with Module 1 

15 using the coupled equilibria calculations. Note the close agreement between 
experiments and predictions in the bottom table and the disagreement between 
experiments and the predictions using the naive simple hybridization calculation (top 
table of Figure 8). Also note the good agreement in the bottom table for the fully 
matched A-T sequence and mismatch A- A, A-C, and A-6 sequences, thus validating 

20 the mismatch parameters. 

Further, the net hybridization calculations can be extended to 
different temperatures as shown in Figure 9, to reveal how the concentrations of all 
species change with temperature. Given the extinction coefficients and fluorescence 
quantum yields, the concentration vs. temperature profiles shown in Figure 9 can 
25 be used to calculate the fluorescence vs. temperature profile (not shown), thereby 
allowing the prediction of the temperature which produces the maximum 
fluorescence signal and minimum background fluorescence signal. 
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Another manifestation of the concentration calculations is for match 
vs. mismatch discrimination (Figure 10), whereby the concentrations of all species 
at all temperatures can be calculated (Figure 11). For the particular case shown, 
optimal match vs. mismatch discrimination is predicted to occur at 0°C. The 
5 concentration calculations can be generalized for cases in which molecules can form 
many different competing unimolecular, biomolecular, and higher order complexes 
(Figure 12) using generalized equations such as shown in Figure 13 for the five 
molecule case, and solved using the algorithm in Figure 7. 

Algorithm 

10 The hybridization prediction algorithm of the present invention is 

based on a nearest-neighbor-model analysis of the sequences. The algorithm 
accounts for structural motifs including Watson-Crick base pairs (Allawi and 
SantaLucia, 1997; SantaLucia, 1998; Sugimoto et al., 1995; Xia et al., 1998), 
single internal mismatches (Allawi and SantaLucia, 1997; Allawi and SantaLucia, 

15 1998; Allawi and SantaLucia, 1998; Allawi and SantaLucia, 1998; Kierzek et al., 
1999; Peyret et al., 1999; SantaLucia, 1998), double mismatches (Allawi and 
SantaLucia, 1997) coaxial-stacking interfaces (disclosed herein) (Walter and Turner, 
1994), terminal mismatches (disclosed herein) (Freier et al., 1986) and dangling 
ends (Bommarito et al., 2000; Freier et al., 1986). Once the motifs are identified 

20 and their thermodynamic contributions are added, the sum may be corrected for salt 
effects (sodium and magnesium) and the net hybridization is calculated when 
appropriate. 

Algorithm Functions 

A first or main module of the algorithm calculates the hybridization 
25 thermodynamics (AH°, AS°, AG° 37 , T m ) of a given duplex. Net hybridization 
accounting for secondary structure in both strands is also calculated. 
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Parameterization 

Parameters are organized in three arrays. The first array contains 
internal element parameters: Watson-Crick nearest neighbors and single mismatch 
nearest neighbors. The second array contains terminal element parameters: terminal 
5 mismatches and dangling ends. A single parameter is used to account for double 
mismatches except for tandem G T mismatches, which are explicitly enumerated 
(Allawi & SantaLucia, 1997). The third array contains coaxial-stacking parameters 
(contained herein). 

For DNA sequences, the thermodynamic contribution of all Watson- 
10 Crick nearest neighbors and single internal mismatches has been systematically 
studied (Allawi and SantaLucia, 1997; Allawi and SantaLucia, 1998; Allawi and 
SantaLucia, 1998; Allawi and SantaLucia, 1998; Peyret et al., 1999). A limited 
number of sequences containing double mismatches has also been studied (Allawi 
and SantaLucia, 1997). The contributions of dangling ends (Bommarito et al., 
15 2000) have also been systematically analyzed. Salt corrections are available for 
sodium in the range 0.01 to 1 M (SantaLucia, 1998). 

For RNA sequences, the thermodynamic contribution of all Watson- 
Crick nearest neighbors has been systematically studied (Xia et al., 1998). A limited 
number of sequences containing single mismatches has also been studied (Kierzek 
20 et al., 1999). The contribution of dangling ends and terminal mismatches has also 
been systematically analyzed (Freier et al., 1986). No salt correction has been 
developed for RNA and therefore the DNA salt corrections are assumed. These 
corrections are likely to be deficient in the case of RNA. 

For DNA/RNA hybrids, the thermodynamic contribution of all 
25 Watson-Crick nearest neighbors has been systematically studied as well as a limited 
number of sequences containing single mismatches (Sugimoto et al., 1995). As no 
salt correction has been developed for DNA/RNA hybrids, the DNA corrections are 
assumed. The applicability of these corrections to DNA/RNA hybrids has not been 
tested. The parameter arrays are designed to easily accommodate implementation 
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of new parameters and salt corrections including thermodynamics parameters for 
modified bases and denaturant effects. 

Correction for Hybridization to DNA Microchips 

A linear correction of the free energy is implemented in the algorithm 
5 of the invention to correct for hybridization to DNA microchips: 

AG° 37 (microchip) = aAG° 37 (solution) +b (32) 

where a and b are user defined real coefficients. Fotin et aL (Fotin et al., 1998) 
showed that a linear relationship could be used to relate the free energies obtained 
for hybridization in solution and on microchip surfaces. However, the relation 
10 between thermodynamics measured in solution and thermodynamics measured using 
microarrays is still unclear and appears to be different depending on the 
manufacture and type of microarrays. 

User Interface: Input and Output 

Figure 2a shows the user interface input. The users enter the 
15 sequence of each strand, the hybridization conditions (hybridization temperature, 
strand concentrations, and monovalent cations and concentrations), and 
thermodynamic corrections for single strand folding. Figure 2b shows the output 
corresponding to the input in Figure 2a. 

The algorithm can be used via the Internet at: 
20 http://isll.chem.wavne.edu/Hvther/hvthermlmain.html . The algorithm may be 
written in FORTRAN 77 and run on UNIX environment or other languages and 
environments. 
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Molecular Beacons 

The algorithm may be used to predict the thermodynamics of a set 
of literature measurements for molecular beacons (Bonnet et al., 1999). Molecular 
beacons are high specificity probes that are efficient for mutation analysis 
5 (Giensendorf et aL, 1998) and multiplex detection of single nucleotide variations 
(Marras et al., 1999). The design and efficiency optimization of these beacons is 
helped by hybridization thermodynamics prediction. Bonnet et al. studied, the 
hybridization of the molecular beacon 5 CGC, TCC, CAA, AAA, AAA, AAA, 
CCG AGC G 3 to a set of four different targets including a perfect match duplex, 
10 and three different duplexes containing one mismatch. Free energy and enthalpy for 
duplex folding may be calculated using the DNA MFOLD program 
(http://mfold2.wustl.edu/-mfold/dna/forml.cgi). These parameters may then 
incorporated as secondary structure corrections in Figure 2a. 

The software to implement the algorithm may be written in 
15 FORTRAN, C + + , Visual Basic, HTML, and JAVA script computer languages. 
Two graphical user interfaces may be provided: Windows application and web 
browser format. The software may run on IBM/PC, Sun, and Silicon Graphics 
platforms. 

The software may be written in several modules as described below. 

20 A. Interactive Mode: Command Line Interface in MS-DOS 

MODULE 1 (As Previously Described above) 

Function . Module 1 predicts the hybridization thermodynamics of a given 
duplex (DNA/DNA, RNA/RNA, or DNA/RNA). 
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Input (Figure 2a) 

Input of Sequences 

1. Only the following characters are accepted: A, a, C, c, G, g, 
T, t, U, u, /, *, +. Single blank characters and numbers will be automatically 

5 edited, but more than one carriage return is not permitted. 

2. If the duplex contains a dangling end on a strand, the 
sequence of the other strand should contain a * at the corresponding position. (This 
is very important to include for primer binding to a large target sequence). Note: 
The top strand must be entered in 5 1 to 3' orientation, but the bottom strand must 

10 be entered in 3' to 5 1 orientation. Also, a " + " must be added at the end of each 
sequence. There is a length limit of 1024 characters for sequence entries. In 
module 1, it is important to be sure that both sequences have the same length. 

Example: AAAACCCCTGA + 
*TTTGGGGAC* + 

15 3. Only the bottom strand may contain coaxially stacked 

nucleotides. A "/" should be inserted at the site of a strand nick (i.e. between the 
coaxially stacked nucleotides). This feature is useful for predicting stacked 
hybridization stability. 

Example: AAAACCCCC + 
20 TTTT/GGGG+ 

Input of Salt and Strand Concentrations 

The monovalent salt should be the sum of all monovalent cation 
concentrations in a solution in units of molarity. For example, a solution of 100 
mM KC1, 50 mM NaCl, 10 mM Na 2 P0 4 , 0.1 mM Na 2 EDTA would account for a 
25 total of 0. 1702 M monovalent. The thermodynamic predictions are applicable over 
a salt range of 0.01 to 1 M monovalent cation. The correction applied is from 
SantaLucia (1998) Proc. Natl Acad. Sci. 95, 1460. The sodium correction applies 
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for oligonucleotides with fewer than about 30 base pairs. For longer duplexes a 
polymer correction is required, but this is not currently implemented. 

Strand concentrations are entered in units of molarity. The program 
will accept virtually any physically relevant strand concentration. 

5 Hybridization temperature is in Celsius degrees. The limits are 0 to 

100 degrees. 

Special corrections for single-stranded secondary structure and for 
surface corrections for hybridization arrays can be input. The units for input AG° 
are kcal/mol. To determine estimates of single-strand folding energies, see Michael 
10 Zuker's RNA or DNA-MFOLD servers (see http: 
//mfold2.wsutl.edu/~mfold/dna/forml .cgi). The current thermodynamic prediction 
software incorporates the special corrections for single-stranded secondary structure 
and for surface corrections for hybridization arrays. 

For DNA chip arrays, a linear correction can be applied. The user 
15 inputs the slope and intercept coefficients. Based on the work of Mirzabekov 
group, a slope of +1.1 and intercept of +3.2 are appropriate (see Fotin et al. 
(1998) Nucleic Acids Res. 26, 1515-1521). 

Output (Figure 2b) 

Module 1 outputs the hybridization thermodynamics at 1.0 M NaCl 
20 and 37°C (the conditions under which the thermodynamic predictions are most 
accurate), under the salt temperature conditions specified by the user, and also 
displays the net hybridization Tm and AG° if the user specifies that special 
corrections are needed (this allows for single-strand secondary structure of both the 
target and probe DNA to be accounted and for surface effects of chip arrays). 
25 Predictions of AG°, AH°, AS°, and Tm are provided. 
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MODULE 2 

Function . Module 2 finds the best primers of given length 
complementary to a long target nucleic acid. DNA/DNA, RNA/RNA, DNA/RNA 
hybridization types are accepted. The user selects the number of primers to output, 
5 and the program finds the most stable primers and gives their hybridization position 
and thermodynamics of each primer. 

Input (Fisure 3a) 

Input of Salt and Strand Concentrations 

The input of strand and salt concentrations is similar to Module 1 . 

10 Input of Sequences 

The target sequence is input as in Module 1 . 

Output (Figure 3b) 

Primer Length and Number of Best Primers 
Module 2 displays "number of best primers" best primers of length 
15 "primer length" in order of decreasing stability. 

Output 

Module 2 outputs "number of best primers" best primers of length 
"primer length" in order of decreasing stability along with their hybridization 
thermodynamics . 

20 MODULE 3 

Function . Module 3 walks a given primer along a given target and 
finds the thermodynamics for the best target/primer complex and for the competitive 
target/primer complexes: DNA/DNA, RNA/RNA, DNA/RNA, hybridization types 
are accepted. 
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Input (Figure 4a) 

Input of Sequences 

The input is similar to Module 1 . The target has to be longer than 

the primer. 



Percent Stability p of Alternative Binding Sites 
Compared to the Most Stable Binding Site 

This parameter excludes all competitive sites that are not within the 
10 defined percent of the best primer stability. If the best primer stability is -5 
kcal/mol and p= 10 then any competitive site of energy higher than -5 + (10/100*5) 
= -4.5 kcal/mol will not be displayed. 

Number of Base Pairs Required to Compute the Solution 
This parameter excludes all competitive sites that contain less 
15 Watson-Crick base pairs than the defined value. 

Output (Figure 4b) 

Module 3 outputs the best primer binding site and the competitive 
binding sites that pass the filtering criteria (percent stability p of alternative binding 
sites compared to the most stable binding site and number of best primers). 

20 MODULE 4 



Input of Salt and Strand Concentrations 

The input of salt and strand concentrations is similar to Module 1 . 



Function , Batch mode calculations (see below). 
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MODULE 5 

Function , Module 5 is a combination of Modules 2 and 3 and finds 
the n best primers of given length complementary to a given section of a target and 
display the thermodynamics of the target/primer system(s). Then, each best primer 
5 is walked along the whole target to find the competitive hybridization sites. The 
thermodynamics of the target/primer systems at these alternative sites is then 
displayed, DNA/DNA, RNA/RNA, DNA/RNA, hybridization types are accepted. 

Input (Figure 5a) 

Input of Sequences 
10 The target sequence is input as in Module 1 . 

Input of Salt and Strand Concentrations 

The input of salt and strand concentrations is similar to Module 1 . 

Sequence Section Where to Find the Best Primers 
Module 5 finds the best primers in the target region ranking from 
15 "position of initial nucleotide" to "position of final nucleotide" . Note that Module 
5 then looks for competitive sites of each best primers in the whole target. 

Percent Stability of Alternative Binding Sites 
Compared to the Most Stable Binding Site 

The function of this parameter is the same as in Module 3. This 
20 parameter is input for each best primer corresponding to the "number of best 
primer" specified. 

Number of Base Pairs Required to Compute the Solution 
The function of this parameter is the same as in Module 3. This 
parameter is input for each best primer corresponding to the "number of best 
25 primer" specified. 
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Output (Figure 5b) 

Primer Length and Number of Best Primers 
Module 5 displays "number of best primers" best primers of length 
"primer length" by order of decreasing stability. 

5 Output 

Module 5 displays "number of best primers" best primers and their 
competitive sites by order of stability along with their hybridization 
thermodynamics. The best primer and its ranked competitive hybridization sites are 
listed first. Then, the second best primer is listed with its competitive hybridization 
10 sites. 

MODULE 6 

Function . Module 6 is similar to Module 3 and walks a given primer 
along a given target and finds the thermodynamics for the best target/primer 
complex and for the competitive target/primer complexes: DNA/DNA, RNA/RNA, 
15 DNA/RNA, hybridization types are accepted. Then, Module 6 simulates the 
concentration of every species at every degree from 1 to 100 °C, as illustrated in 
Figure 6. 

Input (Not Shown) 

Input of Sequences 

20 The input is similar to Module L The target has to be longer than 

the primer. 

Input of Salt and Strand Concentrations 

The input of salt and strand concentrations is similar to Module 1 , 
Percent Stability p of Alternative Binding Sites 
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Com pared to the Most Stable Binding Site 

This parameter excludes all competitive sites that are not within the 
defined percent of the best primer stability. If the best primer stability is -5 
kcal/mol and p=10, then any competitive site of energy higher than: -5 + 
(10/100*5) = -4.5 kcal/mol will not be displayed. 

Number of Base Pairs Required to Comput e the Solution 
This parameter excludes all competitive sites that contain less 
Watson-Crick base pairs than the defined value. 

Correction for Target/Target Interaction, Target folding, 
Primer/Primer Interaction and Primer Folding 

The user is asked if he wants to correct for the interactions above. 

If the answer is "y", the user is prompted for AH° 37 corresponding to the 

interaction. Secondary structure thermodynamics can be determined using the 

Zuker algorithm as discussed in Module 1 section. 

Output (Not Shown) 

Concentration Output Filename 

The results from the concentration simulations (concentration of 
species at every temperature) are saved in this file. 

Output 

Module 6 outputs the best primer binding site and the competitive 
binding sites that pass the filtering criteria (percent stability p of alternative binding 
sites compared to the most stable binding site and number of best primers). The 
concentration simulations are saved in a file specified by the user. 
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MODULE 7 

Function . Module 7 is a combination of Modules 2 and 5 and finds 
the n best primers of given length complementary to a given section of a target and 
display the thermodynamics of the target/primer system(s). Then, each best primer 
5 is walked along the whole target to find the competitive hybridization sites. The 
thermodynamics of the target/primer systems at these alternative sites is then 
displayed. DNA/DNA, RNA/RNA, DNA/RNA hybridization types are accepted. 
Then, Module 7, like Module 6, simulates the concentration of every species at 
every degree from 1 to 100°C, as illustrated in Figure 6. 

10 Input (Not Shown) 

Input of Sequences 

The target sequence is input as in Module 1 . 
Input of Salt and Strand Concentrations 

The input of salt and strand concentrations is similar to Module 1 . 

15 Sequence Section Where to Find the Best Primers 

Module 7 finds best primers in the target region ranking from 
"position of initial nucleotide" to "position of final nucleotide. " Note that Module 
7 then looks for competitive sites of each best primers in the whole target. 

Percent Stability of Alternative Binding Sites 
20 Compared to the Most Stable Binding Site 

The function of this parameter is the same as in Module 3 . This 

parameter is input for each best primer corresponding to the "number of best 

primer" specified. 
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Number of Base Pairs Required to Compute the Solution 

The function of this parameter is the same as in Module 3 . This 

parameter is input for each best primer corresponding to the "number of best 

primer" specified. 

5 Correction for Target/Target Interaction, Target folding, 
Primer/Primer Interaction and Primer Folding 

For each best primer, the user is asked if he wants to correct for the 

interactions above. If the answer is "y", the user is prompted for AH° and AG° 37 

corresponding to the interaction. Secondary structure thermodynamics can be 

10 determined using the Zuker algorithm as discussed in Module 1 section. 

Concentration Output Filenames 

For each best primer, the results from the concentration simulations 
(concentration of species at every temperature) are saved in this file. The user has 
to select a different filename for each best primer. 

15 Output (Not Shown) 

Output 

Primer Length and Number of Best Primers 

Module 7 displays "number of best primers" best primers of length 
"primer length" by order of decreasing stability. 

20 Module 7 displays "number of best primers" best primers and their 

competitive sites by order of stability along with their hybridization 
thermodynamics. The best primer and its ranked competitive hybridization sites are 
listed first. Then, the second best primer is listed with its competitive hybridization 
sites. For each best primer, a file named by the user contains the concentration 

25 simulations. 
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Module 7 allows the user to design optimal primers for applications 
where multiple simultaneous hybridization reactions are occurring, including match 
vs. mismatch hybridization, molecular beacons, DNA oligonucleotide arrays, and 
multiplex PCR. 

5 One commercially important example for the use of Module 7 for 

primer design in a complex hybridization solution is Multiplex PCR, as shown in 
Figure 7. Module 7 allows the user to design optimal primers for Multiplex PCR 
where multiple primers have equal stabilities in binding to the target DNA. Several 
primers must be designed to specifically bind to different sites on target DNA at a 
10 given temperature with minimal background binding to mismatch sites and with 
minimal cross-hybridization between pairs of primers. 

Module 7 minimizes potential primer dimer formation and mismatch 
hybridization for all combinations of input primers. Module 7 optimizes primer 
sequence position, length, and concentration for each primer in relation to all other 
15 species in solution and provides a hybridization profile at all temperatures from 0 
to 100 °C. 

Batch Mode 

MODULE 4 

Function . Module 4 allows any of the previous modules to be run 
20 in batch mode using text files to submit the input and having the data output as text 
files also. 

Type of Input Files 

There are two types of input files: 1) parameter input file, and 2) 
sequence input file. Parameter input files describe what modules to run with what 
25 hybridization parameters and on how many sequences to run them. Example of 
parameter input files for each module with comments are given in the "Batch mode 
parameter files folder." Sequence files contain the sequences that are going to be 
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hybridized in the conditions described by the parameter input files. Examples of 
parameter input files for each module with comments are given in the "Batch mode 
sequence files folder." 

Note that a parameter file can successively run different modules on 
5 various different sequences. 

The user is successively asked for the names of the parameter input 
file, the sequence input file and the thermodynamic data output file. Note that these 
files have to be in the directory containing the executable version of the software. 
Output files will also be created in this same directory. Names of the concentration 
10 simulation files are specified in the parameter input files. 

Examples of Batch Mode Parameter Files 

Comments in parentheses describe the meaning of each entry (note 
that an actual parameter file must not contain these comments). 

DUP (Module 1 : Simple duplex calculations) 

15 1 (Number of sequences to apply this parameter file to) 

1 (Monovalent cations concentration mol/L)) 

1 Mg 2+ concentration mol/L) 

37.0 (Hybridization temperature) 

4e-4 (Top strand concentration mol/L) 

20 4e-4 (Bottom strand concentration mol/L) 

1 (Correction for microchips: slope) 

0 (Correction for microchips: intercept) 

0 (Correction for top strand folding: AG° 37 ) 

0 (Correction for top strand folding: AH° 37 ) 

25 0 (Correction for bottom strand folding: AG° 37 ) 

0 (Correction for bottom strand folding: AH° 37 ) 

END (End of file required) 
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NBP (Module 2: N-best primers) 

1 (Number of sequences to apply this parameter file to) 

1 (Monovalent cations concentration mol/L) 

1 (Mg 2+ concentration mol/L) 

5 37 (Hybridization temperature) 

4e-4 (Top strand concentration mol/L) 

4e-4 (Bottom strand concentration mol/L) 

4 (Primer Length) 

3 (Number of best primers) 

10 1 (Correction for microchips: slope) 

0 (Correction for microchips: intercept) 
END (End of file required) 

PWA (Module 3 primer walk match vs. mismatch sites identification) 

1 (Number of sequences to apply this parameter file to) 
15 1 (Monovalent cations concentration mol/L) 

1 (Mg 2+ concentration mol/L) 
37 (Hybridization temperature) 
4e-4 (Top strand concentration mol/L) 
4e-4 (Bottom strand concentration mol/L) 

20 90 (Percent window of best primer stability for alternative sites) 

2 (Number of WC base pairs required to compute the solution) 
cgcg + (Primer sequence, + required) 

1 (Correction for microchips: slope) 

0 (Correction for microchips: intercept) 
25 END (End of file required) 

BPW (Module 5 displays "number of best primers" best primers and their 
competitive sites by order of stability along with their hybridization 
thermodynamics) 

1 (Number of sequences to apply this parameter file to) 
30 1 (Lower limit of primer search area) 

10 (Upper limit of primer search area) 
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1 Mg 2+ concentration mol/L) 

37 (Hybridization temperature) 

4e-4 (Top strand concentration mol/L) 

4e-4 (Bottom strand concentration mol/L) 

5 4 (Primer length) 

1 (Number of best primers) 

1 (Correction for microchips: slope) 

0 (Correction for microchips: intercept) 

800 (Percent window of best primer stability for alternative sites) 

10 2 (Number of WC base pairs required to compute the solution) 

END (End of file required) 

PWC (Module 6 primer walk with concentration calculations) 

1 (Number of sequences to apply this parameter file to) 
1 (Monovalent cations concentration mol/L) 

15 1 Mg 2+ concentration mol/L) 

37 (Hybridization temperature) 

4e-4 (Top strand concentration mol/L) 

4e-4 (Bottom strand concentration mol/L) 

90 (Percent window of best primer stability for alternative sites) 

20 2 (Number of WC base pairs required to compute the solution) 

cgcg+ (Primer sequence, + required) 

1 (Correction for microchips: slope) 

0 (Correction for microchips: intercept) 

0 (Correction for target folding: AG° 37 ) 

25 0 (Correction for target folding: AH°) 

0 (Correction for target/target interaction: AG° 37 ) 

0 (Correction for target/target interaction: AH°) 

0 (Correction for primer folding: AG° 37 ) 

0 (Correction for primer folding: AH°) 

30 0 (Correction for primer/primer interaction: AG° 37 ) 

0 (Correction for primer/primer interaction: AH°) 

outconc (Concentration output file name) 

^37- 
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END (End of file required) 





BWC 


(Module 7, N-best primers, primer walk, and concentration 






calculations) 




1 


(Number of sequences to apply this parameter file to) 


5 


1 


(Lower limit of primer search area) 




10 


(Upper limit of primer search area) 




1 


(Monovalent cations concentration mol/L) 




0 


Mg 2+ concentration mol/L) 




37 


(Hybridization temperature) 


10 


4e-4 


(Top strand concentration mol/L) 




4e-4 


(Bottom strand concentration mol/L) 




4 


(Primer length) 




1 


(Number of best primers) 




1 


(Correction for microchips: slope) 


15 


0 


(Correction for microchips: intercept) 




800 


(Percent window of best primer stability for alternative sites) 




2 


(Number of WC base pairs required to compute the solution) 




0 


(Correction for target folding: AG° 37 ) 




0 


(Correction for target folding: AH°) 


20 


0 


(Correction for target/target interaction: AG° 37 ) 




0 


(Correction for target/target interaction: AH°) 




n 
\j 


^v^urrctAiun lor pruner ioiuing. hkj 37 j 




0 


(Correction for primer folding: AH°) 




0 


(Correction for primer/primer interaction: AG° 37 ) 


25 


0 


(Correction for primer/primer interaction: AH°) 




outconc 


(Concentration output file name) 




END 


(End of file required) 




PPW 


(Module 8: walk a primer along itself to find interaction sites. PWA 






with probe = primer) 


30 


1 


(Number of sequences to apply this parameter file to) 




1 


(Monovalent cations concentration mol/L) 
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1 


Mg 2+ concentration mol/L) 


37 


(Hybridization temperature) 


4e-4 


(Primer concentration mol/L) 


900 


(Percent window of best primer stability for alternative sites) 


2 


(Number of WC base pairs required to compute the solution) 


END 


(End of file required) 



Examples of Batch Mode Sequence Files 

For Module 1 : Pup 
1 (Sequence number) 

10 agcgca-h (Top strand sequence) 

tcgcgt+ (Bottom strand sequence) 

For Module 2: NBP 

1 (Sequence number) 

agcgca+ (Target sequence) 

15 For Module 3: 

1 (Sequence number) 

cgcctgcggccc+ (Target sequence) 

For Module 5: bpw 
1 (Sequence number) 

20 cgcctgcgccc-h (Target sequence) 

For Module 6: pwc 

1 (Sequence number) 

agcgcaH- (Target sequence) 

For Module 7: bwc 
25 1 (Sequence number) 

agcgcaH- (Target sequence) 
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For Module 8: ppw 

1 (Sequence number) 

agcgca+ (Primer sequence) 

Example of Batch Mode Parameter and Sequence Files 
5 to Run Different Modules Successively 



Parameter File: 



10 



DUP 
2 

0.05 

1.5e-3 

37.0 

le-6 

2e-7 

1 



15 0. 



-2.12 
-37.3 



(Executes Module 1) 

(Apply to Module 1 to 2 sequence sets) 



0 

20 PWC 
1 

0.16 

0.0025 

37 

25 10e-9 
le-9 
800 
8 

TCGAACGTAC+ 
30 1 



(Executes Module 6) 

(Apply to Module 6 to 1 sequence set) 
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0 
0 
0 
0 

5 0 
0 
0 
0 
0 

10 outwash 

DUP (Executes Module 1) 

4 (Apply to Module 1 to 4 sequence sets) 

1 

0 

15 37.0 
le-6 
le-6 
1 
0 

20 0 
0 
0 
0 

END 

25 Other modules can be similarly appended. 
Sequence File 

1 (input for Module 1) 

ttgcctaggggaccaggtccaact + 
aacggatcccctggtccaggttga + 
30 2 (input for Module 1) 
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ttgcctaggggaccaggtccaact + 
aacggatcccctggtccaggttga + 
3 

C AGCTTGC ATG A A A AGCTTGCGTGT + 
5 4 

AAAAAA+ 

TTTTTT+ 

5 

acgcgc + 
10 tgcgcg + 
6 

gggaaagggg + 
*cctttccc* + 
7 

15 tttaaattt-h 
aaatttaaa+ 
8 

cgcgtgagggcc+ 
gcgctctccccgg+ 



(input for Module 6) 
(input for Module 1) 

(input for Module 1) 

(input for Module 1) 

(input for Module 1) 

(input for Module 1) 



20 Parameterization of the Algorithm of the Invention 

Caution: RNA/RNA and DNA/RNA duplexes contain motifs for 
which no literature data are available. In these cases, DNA/DNA parameters are 
assumed. Therefore, predictions might be inaccurate. Users are encouraged to use 
this program with caution and discernment. 
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No data are available for the following motifs: 



RNA/RNA single mismatches 
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RNA/DNA single mismatches 
dangling ends 
terminal mismatches 
Single mismatches 

5 Double mismatch parameters are estimated for all types of duplexes 

(DNA/DNA, RNA/RNA, DNA/RNA). 





DNA THERMODYNAMIC PARAMETERS 




Watson-Crick nearest-neighbors 






10 


SantaLucia, Allawi, and Seneviratne 
(1996) Biochemistry 35, 3555; 
Allawi and SantaLucia (1997) 
Biochemistry 36, 10581 


12 parameters 


108 sequences 




Single mismatch nearest-neighbors 






15 


Allawi and SantaLucia (1997) 
Biochemistry 36, 100581; 
Allawi and SantaLucia (1997) 
Nucleic Acids Res., 26, 2694; Peyret 
et al., (1999) Biochemistry 38, 3468 


44 parameters 180 sequences 

Allawi and SantaLucia (1998) Biochemistry 37, 2170 
Allawi and SantaLucia (1998) Biochemistry 37, 9435 


20 


Terminal Mismatch nearest- 
neighbors 

(Appendix) 


48 parameters 


48 sequences 




Dangling end nearest-neighbors 






25 


S. Bonunarito, Peyret, SantaLucia 
(2000) Nucleic Acids Res. 28, 1929- 
1934. 


16 parameters 


16 sequences 




Na + dependence law 








SantaLucia (1998) Proc. Natl Acad. 
Sci. USA 95, 1460-1465 


1 parameter 


86 sequences 
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RNA THERMODYNAMIC PARAMETERS 


Watson-Crick nearest-neighbors 

Xia et at. (1998) Biochemistry 37, 14719 


12 parameters 


Sin<r1f* \f inmate* Tips 

Kierzek, Burkard, and Turner (1999) 
Biochemistry 38, 14214-14223 


44 parameters 


XwlilllllCLl llllOllluLwIi ilbalwQl 11V 1 U Ul o 

Freier et al. (1986) Proc. Natl Acad. Sci. USA 83, 
9373 


48 parameters 


Dangling end nearest-neighbors 

Freier et al. (1986) Proc. Natl Acad. Sci. USA 83, 
9373 


32 parameters 


Coaxial stacking nearest-neiehbors 

Walter and Turner (1994) Biochemistry 33, 12715 


16 parameters 


Loop parameters 

Matthews et al. (1999) /. Mol. Biol 288, 911-940 
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HYBRID DNA/RNA THERMODYNAMIC PARAMETERS 


Watson-Crick nearest-neighbors 

Sugimoto et al. (1995) Biochemistry 34, 11211 


17 parameters 


rU»dG and rG^dT mismatches 

Sugimoto et al. (1997) Nucleic Acids Symp. Ser. 37, 199 
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DNA LOOP THERMODYNAMIC PARAMETERS 



Hairpins 

Hilbers et al. (1985) Biochie 67, 685-695 

Blommers et al. (1989) Biochemistry 28, 7491-7498 

Antaoetal. (1991) Nucleic Acids Res. 19, 5901-5905 

Antao et al. (1991) Nucleic Acids Res, 20, 819-824 

Senior et al. (1988) Proc. Natl. Acad. Sci. USA 85, 6242-6246 



10 



Bulges 

LeBlanc and Morden (1991) Biochemistry 30, 4042-4047 

Zieba et al. (1991) Biochemistry 30, 8018-8026 

Ke et al. (1995) Biochemistry 34, 4593-4600 

Turner, D.H. (1992) Curr. Opin. Struc. Biol 2, 334-337 



15 



Multibranched Loops 

Kadrmas et al. (1995) Nucleic Acids Res. 23, 2122 
Lilley and Hallam (1984) /. MoL Biol 180, 179-200 
Lu et al. (1991) /. MoL Biol. 223, 781-789 
Ladbury et al. (1994) Biochemistry 33, 6828-6833 
Leontis et al. (1991) Nucleic Acids Res. 19, 759-766 



20 
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The parameters for multibranched loops are from a best fit analysis of 
secondary structure predictions vs. experiments as done by Jaeger et al. for 
RNA (Jaeger et al. (1989) PNAS 86, 7706-7710). The current parameters for 
multibranched loops neglect the sequence and complicated length dependence 
described by Leontis and coworkers, but approximate 4-way junctions fairly 
well. Implementation of more complicated rules will require modification of 
the MFOLD algorithm. 



While the best mode for carrying out the invention has been 
described in detail, those familiar with the art to which this invention relates will 
recognize various alternative designs and embodiments for practicing the invention 
as defined by the following claims. 
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Table 1 : Thermodynamic Parameters for Duplex Formation in 1M NaCl a . 











AH° b 


AS ob 


AGV " 












(kcal / mol) 


(cal / mol K) 


(kcal / mol) 


(°C) 


A T 


G A G C 


T C A 


A 


-56.7 ± 2.5 


-155.8 ± 3.7 


-9.07 ±0.12 


52.3 


A A 


P T C G 


A G T 


A 










A A 


G A G C 


T C T 


A 


-55.1 ± 1.4 


-149.2 ± 3.8 


-8.91 ± 0.12 


56.0 




c T r g 

W X w VJ 


A G A 

Ji VJ xm. 


A 










A G 


T A G C 


T A C 


A 


-60.3 ± 1.8 


-165.8 ± 5.1 


-9.03 ±0.12 


54.5 


A C 


A T P A 

jTv X V-' 


A T G 

-TV X VJ 


A 










A C 


G A T A 


T C G 


A 


-67.8 ± 1.3 


-192.0 ± 3.5 


-8.87 ± 0.06 


49.4 


a n 

J± v* 


C T A T 


A G C 

xV VJ V— - 


A 










£ T 


G A G C 


T C A 


£ 


-50.6 ± 1.2 


-136.9 ± 2.8 


-8.15 ± 0.11 


52.9 


C A 


r t r n 


A O T 

^Tl VJ A 


c 










£ A 


G A G C 


T C T 


£ 


-51.4 ± 1.3 


-137.5 ± 3.2 


-8.34 ± 0.14 


56.9 


C T 


c t c G 

Vv 1 Vv vj 


AHA 

xv VJ irV 












£ G 


T A G C 


T A C 


£ 


-55.8 ± 1.8 


-154.6 ± 4.8 


-8.04 ± 0.16 


49.7 




A T r G 

x\- X \s VJ 


a t a 












£ C 


G A T A 


T C G 


£ 


-59.8 ± 1.2 


-166.9 ± 3.0 


-8.12 ± 0.07 


49.7 




r T A T 


a n c 

A U V- 












GT 


G A G C 


T C A 


£ 


-52.6 ± 1.3 


-140.7 ± 3.3 


-8.57 ± 0.13 


57.7 


n a 


r t r g 

\-r A V^ VJ 


a n t 

n VJ X 


Li 










£ A 


G A G C 


T C T 


a 


-52.3 ± 1.4 


-142.5 ± 4.1 


-8.34 ± 0.08 


52.3 


G T 


C T C G 


A G A 


G 










£G 


T A G C 


T A C 


£ 


-59.2 ± 1.0 


-164.2 ± 2.8 


-8.68 ± 0.07 


51.2 


GC 


A T C G 


A T G 


£ 










£C 


G A T A 


T C G 


£ 


-65.7 ± 0.9 


-183.6 ± 2.2 


-8.81 ± 0.07 


52.2 


£G 


C T A T 


A G C 


£ 
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Table 1: Continued 3 . 











AH ob 


AS ob 


AG V b 


r-r- C 

T M 










flcra] f vc\cX\ 


(cal / mol K) 


(kcal / mol) 


(°C) 


T T 


G A G C 


T C A 


I 


-55.4 ± 1.1 


-149.4 ± 3.0 


-8.62 ±0.10 


56.9 


I A 


C T C G 


A G T 


I 










I A 


G A G C 


T C T 


I 


-56.5 ± 1.3 


-154.3 ± 3.7 


-8.72 ± 0.09 


54.3 


I T 


C T C G 


A G A 


1 










1 G 


T A G C 


T A C 


1 


-63.8 ± 0.8 


-178.0 ± 2.1 


-8.75 ± 0.06 


52.1 


1 C 


A T C G 


A T G 


1 










T C 


G A T A 


T C G 


I 


-66.8 ± 0.6 


-187.9 ± 1.6 


-8.60 ± 0.04 


50.5 


1 G 


C T A T 


A G C 


1 










C T 


G A G C 


T C A 


A 


-53.6 ± 1.3 


-145.8 ± 4.0 


-8 42 ± 0 06 


53.6 


A A 


C T C G 


A G T 


£ 










A T 


G A G C 


T C A 


c 


-54.0 ± 1.3 


-144.4 ± 3.2 


-8 92 ± 0 15 


58.8 


£ A 


C T C G 


A G T 


A 










£ G 


T A G C 


T A C 


A 


-56.8 ± 1.4 


-155.6 ± 3.5 


-8 53 ± 0 13 


53.7 


AC 


A T C G 


A T G 


Q 










A G 


T A G C 


T A C 


c 


-57.1 ±1.2 


-1S6 0 ± 'SO 




S4 A 


c c 


A T C G 


A T G 


A 










£ C 


G A T A 


T C G 


A 


-61 8 ± 0 5 


-172.7 ± 1.4 


-8 30 ± 0 03 


50.6 


A G 


C T A T 


A G C 


£ 










A C 


G A T A 


T C G 


£ 


-58.3 ± 1.7 


-158.5 ± 4.2 


-8.91 ± 0.17 


56.8 


£ G 


C T A T 


A G C 


A 










£ A 


G A G C 


T C T 


A 


-54.6 ± 0.6 


-147.9 ± 1.4 


-8.66 ± 0.07 


55.4 


AT 


C T C G 


A G A 


£ 
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Table 1 : Continued 3 . 



















AH ob 


AS ob 


AG° 3 , b 






















\y<Xi i iriui 




(°C) 


A A 


G A 


G 


c 


T 


c 


T 


c 


-55.5 ± 1.4 


-150.4 ± 4.4 


-8.85 ± 0.08 


55.8 


£ T 


C T 


C 


G 


A 


G 


A 


A 










T T 


G A 


G 


c 


T 


c 


A 


c 


-52.2 ± 1.0 


-140.8 ± 2.4 


-8.38 ±0.11 


54.9 


£ A 


C T 


C 


G 


A 


G 


T 


I 










£ T 


G A 


G 


C 


T 


C 


A 


I 


-55.1 ± 0.8 


-150.5 ± 1.9 


-8.42 ± 0.08 


53.2 


X A 


C 1 




G 


A 


G 


X 


L 










I G 


T A 


G 


C 


T 


A 


C 


£ 


-58.0 ± 1.4 


-159.6 ± 3.7 


-8.39 ± 0.10 


52.9 


£C 


A T 


C 


G 


A 


T 


G 


I 










£ G 


T A 


G 


C 


T 


A 


C 


I 


-59.4 ± 0.9 


-164.6 ± 2.1 


-8.21 ± 0.06 


51.7 


T C 


A T 


C 


G 


A 


T 


G 


£ 










I C 


G A 


T 


A 


T 


C 


G 


£ 


-61.7 ± 1.1 


-170.6 ± 3.2 


-8.33 ± 0.12 


53.7 


£G 


C T 


A 


T 


A 


G 


C 


I 










£ C 


G A 


T 


A 


T 


C 


G 


I 


-57.9 ± 1.1 


-159.5 ± 2.7 


-8.11 ± 0.12 


52.7 


I G 


C T 


A 


T 


A 


G 


C 


£ 










1 A 


G A 


G 


C 


T 


C 


T 


£ 


-55.0 ± 1.4 


-148.0 ± 3.5 


-8.80 ± 0.12 


57.5 


£ T 


C T 


C 


G 


A 


G 


A 


1 










£ A 


G A 


G 


C 


T 


C 


T 


I 


-51.5 ± 1.1 


-137.9 ± 2.6 


-8.44 ± 0.15 


56.4 


I T 


C T 


C 


G 


A 


G 


A 


£ 










£T 


G A 


G 


C 


T 


C 


A 


A 


-54.2 ± 1.4 


-146.1 ± 3.3 


-8.77 ± 0.16 


56.8 


A A 


C G 


C 


G 


A 


G 


T 


£ 










A T 


G A 


G 


C 


T 


C 


A 


£ 


-55.4 ± 1.2 


-148.6 ± 3.1 


-9.03 ±0.14 


59.1 


£ A 


C T 


C 


G 


A 


G 


T 


A 
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Table 1 : Continued 3 . 











AH ob 


AS ob 


AG° 37 b 


T„ c 










(kcal / mol) 


(cal / moi K) 


(kcal / mol) 




(jr G 


1 A O C 


T A C* 


A 


Oy,H a: 1 .0 




-S 76 ± 0 1 5 


53. S 


A C 


A T C G 


A T G 


Q 










A G 


T A G C 


T A C 


G 


-63.7 ±1.6 


-174.y ± 4.1 


Q /1£ j. A 12 

-y.4o ± u.io 




QC 


A T C G 


A T G 


A 










G C 


G A T A 


T C G 


A 


/ca y< 1 A o 

-60.4 ± u.o 


-166.6 ± z.l 


O £A _t_ A AQ 

-5.5U ± U.KJy 


JO. J 


A G 


C T A T 


A G C 


G 










A 

A C 


VJ A 1 A 


t c 


Si 


£11-1-1'* 
-Ol.i i 1.-? 


1 £7 < -u 7 /* 
•lO / J ± 




JJ.O 




C T A T 


A G C 


A 










G A 


u A u t 


t 1 t 


A 


ZA A J. 1 1 

-54.U x 1.1 


1 A A Q _u "> "7 


C CO _i_ A 1 A 

-o.o^ ± U. 14 


5 / .y 


AT 


C T C G 


A G A 


G 










A A 


u A O I 


T Z^ 1 T* 
1 C 1 


G 


-54.5 ±1.7 


-14o.U ± 5.0 


o nn i A AA 

-o.yo ± o.oy 


56.3 


£T 


C T C G 


A G A 


A 










T* HP 

1 T 


G A G C 


T C A 


G 


-56.8 ± 0.7 


-155.5 ± 1.7 


-8.64 ± 0.06 


53.5 


G A 


C T C G 


A G T 


I 










G T 


G A G C 


T C A 


I 


-57.4 ± 0.7 


-156.7 ± 2.0 


-8.80 ± 0.05 


54.7 


1 A 


C T C G 


A G T 


G 










X G 


1 A (j C 


i A C 


G 


-59.2 ±1.3 


-161.3 ± 3.6 


-8.93 ± 0.08 


56.2 


GC 


A T C G 


A T G 


I 










G G 


1 A G C 


1 A C 


np 

1 


-64.8 ± 2.6 


-182.3 ± 7.2 


-8.63 ± 0.17 


50.0 


I C 


A T C G 


A T G 


G 










I c 


G A T A 


T C G 


G 


-63.3 ± 0.6 


-176.9 ± 1.6 


-8.42 ± 0.05 


51.2 


GC 


C T A T 


A G C 


1 










GC 


G A T A 


T C G 


I 


-63.7 ± 0.9 


-177.4 ± 2.2 


-8.73 ± 0.12 


52.6 


I G 


C T A T 


A G C 


G 
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Table 1 : Continued 3 . 




AH ob 
(kcal / mol) 


AS ob 
(cal / mol K) 


AG°„ b 
(kcal / mol) 


T M C 
CO 


1AGAGC TCTfi 
gTCTCGAGAl 


-57.8 ± 0.6 


-157.4 ± 1.7 


-8.95 ± 0.03 


55.6 


SAGAGCTCT1 
ITCTCGAGAfi 


-57.3 ± 1.2 


-155.5 dfc 3.3 


-8.95 ±0.10 


56.4 


Core sequences 

C GATATCG d 

nCTAT AGP 


-51.9 ± 0.6 


-145.3 ± 1.4 


-6.89 ± 0.09 


44.1 


GTAGC TAC d 
CATCGATG 


-51.6 ± 0.6 


-143.7 ± 1.3 


-7.01 ± 0.08 


45.6 


AGAGCTCT 
T CTCGAGA 


-50.0 ± 0.7 


-136.5 ± 1.7 


-7.76 ± 0.06 


50.2 


T GAGCTCA 
ACTCGAGT 


-50.5 ± 0.5 


-137.7 ± 1.3 


-7.73 ± 0.04 


50.4 



a The top strand of each duplex is represented in the 5* to 3 1 orientation and the bottom strand 
is shown in the 3' to 5' direction. Terminal mismatch nearest neighbors are represented in bold. 

Mismatches are underlined. b AH 0 , AS°, and AG° 37 are the error-weighted averages of the 1/T M vs. 
In C T plot and curve fit methods in Table SI. Errors reflect the precision of the data (see text). 

c T M calculated using 10" 4 M total strand concentration. d Data from reference (19). 
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Tab.* 2: Neighbor .hemodynamic of 
Ukc-with-Uksba- .erminai mismatches n 1 M NaO 

I 



Dimer 



AA/TA 
TA/AA 
CA/GA 
GA/CA 

AQ/TC 
TQ/AC 
CC/GQ 
GC/C£ 

Afi/Tfi 
TQ/AQ 
CQ/GQ 
GG./CG. 

AI/TT 
TI/AI 
CI/GI 
GI/CI 



AH° 



AS C 



Sequence 1 (kcal/mol) 



(e.u) 



Terminal A* A Mismatches 
.3.1 ±1.3 -7.8 ±2.0 

-2.5 ± 0.8 " 6 - 3 ± 2A 

.4.3 ±1-0 -10.7 ±2.6 

.8.0 ±0.7 -22.5 ±1-9 

Terminal C'C Mismatches 
.0.1 ±0.6 0.5 ±1-5 

.0.7 ±0.7 -1.3 ±1-8 

.2.1 ± 0.9 -5-1 ± 2 - 5 

.3.9 ± 0.7 -10.6 ± 1-7 

Terminal G«G Mismatches 
-1.1 ±0.7 -2.1 ±1-8 

.1.1 ± 0.8 -2-7 ± 2.2 

.3.8 ± 0.6 -9-5 ± 1-5 

-0.7 ±0.5 -19.2 ±1.3 

Terminal T»T Mismatches 
.2.4 ± 0.6 -6-5 ± 1-6 

.3.2 ±0.7 -8.9 ±2.1 

.6.1 ±0.5 -16.9 ±1-2 

.7.4 ±0.4 -21.2 ±1.1 



AG 



o »> 
37 



(kcal/mol) 

.0.67 ± 0.06 
-0.58 ± 0.07 
.1.01 ± 0.07 
.0.99 ± 0.06 

.0.21 ± 0.06 
-0.29 ± 0.07 
-0.52 ± 0.09 
-0.62 ± 0.06 

.0.42 ± 0.07 
-0.29 ± 0.05 
.0.83 ± 0.05 
-0.96 ± 0.06 

.0.45 ± 0.05 
-0.48 ± 0.05 
.0.87 ± 0.05 
-0.86 ± 0.05 



Gl/L-i. - ■ ~" 

*^c^an^^ Cn ° TS ^ 
calculated from Table 1 using equations 4 and 5 . 
Dimers are given in antiparallel orientation (e.g. AC™ J* 
J * I -:L „nt h t-TA-51 ^matches are underlined. 
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Tab* 3: Nsares.-ne.ghbor <— na™ P—s 
mix ed-bas C tem.in.1 mismartes in 1 M NaCl 



Dimer 
sequence a 



AH° 
(kcal/mol) 



AS 0 
(e.u) 



Terminal A'C Mismatches 



A£/T£ 
AC/TA 
CA'GC 
C£/GA 
GA/CC 

GC/CA 

TAAC 

TQ/AA 

AC/IT 

AI/TC 

CC/GI 

CI/GC 

GC/CI 

GI/CQ 

TQ/AI 

TI/AC 

AA/TQ 

AQ/TA 

CA/GQ 

CQ/GA 

GA/CQ 

GQ/CA 

TA/AQ 

Tfi/AA 



-1.6 ± 0.7 
-1.8 ±0.7 
-2.6 ± 0.8 
-2.7 ± 0.7 
-5.0 ± 0.4 
-3.2 ± 0.9 
-2.3 ± 0.5 
-2.7 ± 0.8 



-4.0 ± 2.1 
-3.8 ± 1-7 
-5.9 ± 1.8 
-6.0 ±1.6 
-13.8 ± 1.0 
-7.1 ±2.2 
-5.9 ± 1.1 
.7.0 ± 2.4 



Terminal OT Mismatches 



-0.9 ± 0.5 
-2.3 ± 0.5 
-3.2 ± 0.8 
-3.9 ± 0.6 
.4.9 ± 0.6 
-3.0 ± 0.6 
-2.5 ± 0.8 
-0.7 ± 0.6 



-1.7 ± 1.4 
-6.3 ± 1.2 
-8.0 ± 2.0 
-10.6 ± 1-2 
-13.5 ± 1.7 
-7.8 ± 1.5 
-6.3 ± 2.0 
-1.2 ± 1.6 



for 

~AG C 37 
(kcal/mol) 



Terminal G« A Mismatches 



-1.9 ± 0.7 
-2.5 ± 0.7 
-3.9 ± 0.8 
-6.0 ± 0.9 
-4.3 ± 0.5 
-4.6 ± 0.7 
-2.0 ± 0.7 
-2.4 ± 0.9 



-4.4 ± 1.8 
-5.9 ± 1.7 
-9.6 ± 2.1 
-15.5 ±2.1 
-11.1 ± 1-3 
-11.4 ± 1.8 
-4.7 ± 1-6 
-5.8 ± 2.7 



-0.35 ± 0.04 
-0.59 ± 0.08 

.0.76 ± 0.07 

,0.85 ± 0.09 

.0.71 ± 0.05 

.1.01 ± 0.10 

.0.45 ± 0.05 
-0.55 ± 0.05 

-0.33 ± 0.06 
-0.35 ± 0.05 
.0.69 ± 0.07 
.0.60 ± 0.05 
-0.72 ± 0.08 
.0.61 ± 0.08 
-0.52 ± 0.07 
.0.34 ± 0.08 

-0.52 ± 0.08 
.0.65 ± 0.07 
-0.88 ± 0.09 
-1.23 ±0.1 
-0.80 ± 0.06 
.1.08 ± 0.09 
.0.53 ± 0.07 
.0.57 ± 0.05 
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Table 3: Continued 


Dimer 


AH° b 


AS ob 


AG° 37 b 


a 

sequence 


(kcal/mol) 


(e.u) 






Terminal G»T Mismatches 




AG/IT 


-J.Z ± U.4 




-0.45 ± 0.04 


AJ/TQ 


-3.5 ± 0.4 


-9.4 ± 1.2 


-0.54 ± 0,03 


C£/GJ 


-3.8 ± 0.7 


-9.0 ± 1.9 


-0.96 ± 0,06 


CJ7GG_ 


-6.6 ± 1.3 


-18.7 ± 3.6 




GQ/CX 


-5.7 ± 0.4 


-15.9 ± 1.0 


-0.76 ± 0.05 


GJ/Cfi 


-5.9 ± 0.5 


-16.1 ± 1.3 


-0.92 ± 0.07 


Tfi/AI 


-3.9 ± 0.5 


-10.5 ± 1.2 


-0.59 ± 0.03 


TT/AG. 


-3.6 ± 0.7 


-9.8 ± 1.9 


-0.59 ± 0.06 



3 Thermodynamic parameters and their corresponding errors 
are calculated from Table 1 using equations 4 and 5. 

b Dimers are given in antiparallel orientation (e.g. A£/TA equals 
5'-AC-3' paired with 3'-TA-5 f ). Mismatches are underlined. 




Table SI: Thermodynamic Parameters for Duplex Formation in 1M NaCl a . 

















3 


AS' 


3 


AG 0 


3^ 


-T- b 

Hi 














(kcal / mol) 


(cal / mol K) 


(kcal / mol) 


(°C) 


A T 


G A G C 


T 


C A 


A 


c 


-56.0 




3.2 


-151.4 


± 


10.1 


-9.07 


-J- 


0.13 


57.0 


A A 


C T C G 


A 


G T 


A 


d 


-57.6 


± 


4.0 


-156.5 


± 


4.0 


-9.11 




0.46 


56.7 


A A 


G A G C 


T 


C T 


A 


c 


-54.6 


± 


1.8 


-146.9 


± 


5.3 


-8.91 


i 


0.12 


57.7 


A T 


C T C G 


A 


G A 


A 


d 


-56.0 


± 


2.3 


-151-6 


± 


5.4 


-8.97 


± 


0.64 


56.4 


A G 


T A G C 


T 


A C 


A 


c 


-59.9 


± 


2.2 


-164.0 


± 


6.6 


-9.03 




0.12 


55.4 


A C 


A T C A 


A 


T G 


A 


d 


-61.3 


± 


3.3 


-168.3 




7.9 


-9.09 


± 


0.83 


55.3 


A C 


G A T A 


T 


C G 


A 


c 


-60.5 


± 


2.4 


-166.5 


± 


7.7 


-8.86 


± 


0.06 


54.3 


AG 


C T A T 


A 


G C 


A 


d 


-70.8 


± 


1.5 


-198.4 


± 


3.9 


-9.28 


± 


0.34 


53.6 


£ T 


G A G C 


T 


C A 


£ 


c 


-52.3 


± 


4.0 


-142.4 


± 


12.7 


-8.16 


± 


0.11 


52.4 


£ A 


C T C G 


A 


G T 


£ 


d 


-50.4 


± 


1.2 


-136.6 


± 


2.8 


-8.06 


± 


0.36 


52.4 


£ A 


G A G C 


T 


C T 


£ 


c 


-55.6 




2.5 


-152.1 


± 


7.7 


-8.36 


± 


0.14 


52.8 


£ T 


C T C G 


A 


G A 


£ 


d 


-49.9 


± 


1.5 


-134.5 


± 


3.5 


-8.14 




0.45 


53.1 


£G 


T A G C 


T 


A C 


£ 


c 


-55.1 


± 


2.4 


-151.8 


± 


7.2 


-8.04 


± 


0.16 


50.9 


£ C 


A T C G 


A 


T G 


£ 


d 


-56.8 


± 


2.7 


-156.9 


± 


6.6 


-8.08 


± 


0.70 


50.7 


£C 


G A T A 


T 


C G 


£ 


c 


-57.3 


± 


2.6 


-158.4 


± 


8.1 


-8.12 


i 


0.07 


50.8 


£ G 


C T A T 


A 


G C 


£ 


d 


-60.5 


± 


1.3 


-168.3 


± 


3.3 


-8.21 


± 


0.31 


50.9 


G_T 


G A G C 


T 


C A 


£ 


c 


-55.2 


i 


2.1 


-150.4 


± 


6.4 


-8.58 


zfc 


0.14 


54.2 


£A 


C T C G 


A 


G T 




d 


-50.9 


± 


1.7 


-137.1 


± 


3.9 


-8.38 


± 


0.48 


54.4 


£A 


G A G C 


T 


C T 


£ 


c 


-56.8 


± 


2.0 


-155.1 


± 


5.6 


-8.72 


± 


0.26 


54.6 


GT 


C T C G 


A 


G A 


£ 


d 


-48.0 


± 


1.9 


-128.1 


± 


5.9 


-8.30 


± 


0.09 


54.9 
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Table SI: Continued. 













AH° 


AS 0 


AG 0 


r 


T b 












(kcal / mol) 


(cal / mol K) 


(kcal / mol) 


CO 




T A G C 


T A C 




c 


-56.9 


± 


1,5 


-155.4 ± 


4.8 


-8.67 


± 


0.07 


54.2 


GC 


A T C G 


A T G 


£ 




-61.1 




1.4 


-168.5 ± 


3.4 


-8.83 


± 


0.35 


53.9 


ec 


G A T A 


T C G 


G 


c 


6? A 


± 


0 A 


-MO 0 ± 
-i i 


7 4 


-8 79 




0 08 


53.3 


QG 


C T A T 


A G C 


Q 


d 


-66.2 


± 


0.9 


-184.7 ± 


2.3 


-8.93 


± 


0.21 


53.0 


T T 


G A G C 


T C A 


r"T< 
I 


c 


-56.8 


± 


1.5 


-155.1 ± 


4.7 


-8.63 


± 


0.10 




T A 


C T C G 


A G T 


1 


d 


-53.7 




1.7 


-145.6 ± 


3.9 


-8.52 


± 


0.46 


54.4 


1 A 


G A G C 


T C T 


1 




-58.7 




1.6 


-160.3 ± 


4.6 


-8.96 


± 


0.19 


55,4 


T T 


C T C G 


A G A 


I 


d 


-52.9 


± 


2.1 


-142.5 ± 


6.4 


-8.66 


± 


0.10 


55.6 


T G 


T A G C 


T A C 


1 


c 


-63.2 


± 


1.1 


-175.5 ± 


3.4 


-8.75 


± 


0.06 


52.9 




A t c n 


ATM 


T 


d 


-64.4 


± 


1.1 


-179.4 ± 


2.6 


-8.80 


± 


0.25 




1 C 


G A T A 


T C G 


i 


c 


-64.6 


± 


1.3 


-180.4 ± 


4.0 


-8.59 


± 


0.05 


51.7 


1 i_r 


P T A T 


A u L 


JL 


d 


-67.4 




0.7 


-189.4 ± 


1.8 


-8.69 


± 


0.16 


MO 


£ T 


G A G C 


T C A 


A 


c 


-56.0 




2.1 


-153.1 ± 


6.2 


-8.58 


± 


0.15 


53.9 


A A 


C T C G 


A G T 


c 


d 


-52.0 


± 


1.7 


-140.6 ± 


5.3 


-8.39 




0.07 


54.1 


A T 


G A G C 


T C A 


£ 


c 


-57.0 


± 


2.4 


-154.9 ± 


7.1 


-8.94 


± 


0.16 


55,8 


£ A 


C T C G 


A G T 


A 


d 


-52.7 


± 


1.6 


-141.7 ± 


3.6 


-8.73 


± 


0.46 


S6 1 


QG 


T A G C 


T A C 


A 


c 


-57.8 


± 


2.9 


-158.8 ± 


8.9 


-8.54 




0.13 


53.2 


A C 


A T C G 


A T G 


£ 


d 


-56.5 


± 


1.6 


-155.0 ± 


3.8 


-8.46 


± 


0.42 


53.0 


A G 


T A G C 


T A C 


£ 


c 


-59.4 




4.0 


-163.2 ± 


12.2 


-8.74 


± 


0.18 


53.9 


£ C 


A T C G 


A T G 


A 


d 


-56.8 


± 


1.3 


-155.5 ± 


3.1 


-8.60 


± 


0.35 


53.8 


£ C 


G A T A 


T C G 


A 


c 


-61.6 


± 


1.1 


-171.8 ± 


3.5 


-8.30 


± 


0.03 


50.8 


A G 


C T A T 


A G C 


£ 


d 


-61.9 


± 


0.6 


-172.9 ± 


1.5 


-8.30 


± 


0.14 


50.7 
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Table SI: Continued 









AS° 


AG° 37 


T M b 






(kcal / mol) 


(cal / mol K) 


(kcal / mol) 


(°C) 


ACGATATCG£ 


C 


-61.8 ± 


3.2 


-170.3 ± 


9.6 


-8.94 ± 0.19 


54.3 


£ G C T A T A G C A 




-57.0 ± 


1.9 


-155.8 ± 


4.6 


-8.71 ± 0.52 


54.4 


£AGAGCTCTA 


c 


-55.7 ± 


1.1 


-151.6 ± 


3.5 


-8.67 ± 0.07 


54.6 


ATCTCGAGA£ 


d 


-542 ± 


0.7 


-147.2 ± 


1.6 


-8.59 ±0.18 


54.6 


AAGAGCTCT£ 


c 


-JO.O ± 


A 1 

4.1 


-isy.s ± 


iZ.4 


-9.07 ± 0.26 


56.1 


£TCTCGAGAA 


d 


-55.1 ± 


1.5 


-149.1 ± 


4.7 


-8.83 ± 0.08 


55.8 


i 1 U A VJ t 1 LAV,, 


c 


-54.1 ± 


1.9 


-147.4 ± 


5.8 


-o.-4U U.IZ 




CACTCGAGTI 


d 


-51.5 ± 


1.1 


-139.4 ± 


2.6 


-8.27 ± 0.32 


53.5 


£ T GAGC TCAI 


c 


-57.3 ± 


4.5 


-157.5 ± 


14.3 


-8.43 ± 0.09 


52.7 


1ACTCG AGT£ 




-55.0 ± 


0.8 


-150.3 ± 


1.9 


-8.37 ± 0.22 


53.0 


TGTAGC T A C £ 


c 


-jo. / 111 






O.J 


-8.39 ± 0.10 


52.0 


£C ATCG ATGI 


d 


-57.4 ± 


1.9 


-158.2 ± 


4.6 


-8.32 ± 0.48 


52.0 


£ G T A G C TACI 


c 


en ^ t 

-59.3 ± 


1.6 


-160.7 ± 


5.1 


-8.21 ± 0.07 


58.1 


1CATCGATGC 


d 


-59.5 ± 


1.2 


-165.4 ± 


2.3 


-8.19 ± 0.28 


50.7 


1CGATATCGC 


c 


-62.8 ± 


1.2 


-175.7 ± 


3.9 


-8.34 ± 0.12 


50.7 


CGCTAT AGCI 


d 


-57.7 ± 


2.4 


-159.6 ± 


5.7 


-8.16 ± 0.59 


51.0 


CCGATATCGi 


c 


-61.1 ± 


1.8 


-170.6 ± 


5.5 


-8.13 ± 0.13 


50.0 


1GCTAT AGC£ 


d 


-56.3 ± 


1.3 


-155.8 ± 


3.1 


-7.98 ± 0.34 


50.2 


1AGAGCTCT£ 


c 


-56.4 ± 


2.0 


-153.5 ± 


6.3 


-8.81 ± 0.12 


55.2 


CTCTCGAGAI 


d 


-53.8 ± 


1.8 


-145.5 ± 


4.3 


-8.69 ± 0.52 


55.4 


£ A G A G C TCTJ 


c 


-54.8 ± 


2.1 


-149.4 ± 


6.4 


-8.47 ±0.17 


53.7 


1TCTCGAGAC 


d 


-50.3 ± 


1.2 


-135.6 ± 


2.8 


-8.28 ± 0.36 


53.9 
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Table SI: Continued. 



fiTGAGCTCAA 
AACGCGAGTfi 

ATGAGCTCAfi 
fiACTCGAGTA 



AH 0 
(kcal / mol) 

-56.6 ± 2.8 
-53.5 ± 1.6 



-57.7 ± 1.9 
-53.8 ± 1.6 



AS 0 
(cal / mol K) 



AG° 37 
(kcal / mol) 



-154.3 ± 8.4 -8.79 ± 0.17 
-144.6 ± 3.6 -8.64 ± 0.44 



-157.0 ± 5.6 
-145.1 ± 3.6 



-9.05 ±0.15 
-8.85 ± 0.45 



(°C) 

55.1 
55.2 

56.3 
56.5 



SGTAGCTACA c -59.2 ± 4.2 
ATCGATGg d -59.4 ± 1.7 



-162.7 ± 12.9 -8.77 ± 0.16 54.1 
-163.3 ± 4.1 -8.72 ± 0.44 53.8 



AGTAGC T AC £ c -63.3 ± 3.4 -173.4 ± 10.2 -9.46 ± 0.20 56.8 
£C ATCGATGA d -63.8 ± 1.8 -175.2 ± 4.5 -9.47 ± 0.45 56.6 



GCGATATCGA c -62.1 ± 1.3 -172.8 ± 3.8 -8.51 ± 0.10 51.8 

AGCTAT AGCfi d -59.1 ± 1.1 -163.6 ± 2.6 -8.39 ± 0.27 51.9 

ACGATATCGfi c -63.1 ± 2.6 -174.1 ± 7.9 -9.06 ± 0.15 54.6 

fiGCT AT AGC A " -60.4 ± 1.6 -166.0 ± 3.7 -8.92 ± 0.40 54.6 

GAGAGCTCTA c -56.4 ± 1.9 -153.5 ± 5.6 -8.84 ± 0.15 55.4 

ATCTCGAGAfi d -52.8 ± 1.3 -142.4 ± 3.1 -8.67 ± 0.38 55.6 

AAGAGC TCTfi c -57.1 ± 2.6 -154.7 ± 7.6 -9.08 ± 0.20 56.7 

fiTCTCGAGAA d -53.1 ± 2.2 -142.8 ± 6.7 -8.86 ± 0.10 56.8 

ITGAGCTCAfi c -55.2 ± 1.5 -150.0 ± 4.5 -8.63 ± 0.07 54.5 

GACTCGAGT1 d -57.2 ± 0.8 -156.4 ± 1.8 -8.71 ± 0.21 54.3 

GTGAGCTCAI c -57.3 ± 0.7 -156.5 ± 2.3 -8.80 ± 0.05 54.9 

1ACTCGAGTG d -57.6 ± 1.6 -157.4 ± 3.7 -8.81 ± 0.42 54.8 

1GTAGCTACG c -59.9 ± 1.8 -164.5 ± 5.5 -8.93 ± 0.08 54.8 

GCATCGATGI d -58.2 ± 2.0 -158.9 ± 4.8 -8.86 ± 0.53 55.0 

GGTAGCTACI c -62.8 ± 3.6 -174.8 ± 11.0 -8.62 ± 0.17 52.3 

1CATCGATGG 6 -67.1 ± 3.8 -187.8 ± 9.4 -8.79 ± 0.85 52.1 
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Table SI: Continued. 







AH 


0 


AS 0 


AG° 37 


-p b 
A M 






(kcal / mol) 


(cal / mol K) 


(kcal / mol) 


(°C) 


ICGATATCGfi 


C 


-58.9 ± 


2.8 


-163.0 ± 


8.9 


-8.40 ± 0.06 


52.0 


CZ n C T AT A (~l C T 
vj L t 1 n i n u 1/ 1 


d 


-63.5 ± 


0.6 


-177.3 ± 


1.6 




SI 6 


fiCGATATCGI 


C 


6 4- 
-Ol.O 


% 7 


• 1 7ft 4 + 




-8.71 ± 0.14 


53.1 


1GCTAT AGCfi 


d 


-63.9 ± 


0.9 


-177.6 ± 


2.3 


-8.77 ± 0.22 


52.8 


JAGAGCTCTG 


c 


-56.4 ± 


1.5 


-153.3 ± 


4.5 


-8,89 ± 0.06 


55.7 


£T CTCGAGA1 




-58 0 ± 


0 6 


-158 1 ± 


1 9 


-8.96 ± 0.03 


55.6 


fiAGAGC TCTI 


c 


-57.7 ± 


1.6 


-157.2 ± 


4.7 


-8.95 ± 0.10 


55.7 


TT CTCG AGA£ 


d 


-56.6 ± 


1.9 


-153.9 ± 


4.6 


-8.90 ± 0.52 


55.7 


Core sequences 
















C GAT ATCG e 


c 


-55.7 ± 


3.9 


-157.1 ± 


12.1 


-6.93 ± 0.12 


44.1 


VJ V^ 1 1 ix. VJ V^ 


d 


-51.8 ± 


0.6 


-145.1 ± 


1.4 




AA O 


GTAGCTAC e 


c 


-55.1 ± 


2.3 


-155.0 ± 


7.0 


-7.04 ± 0.10 


44.9 


CATCGATG 


d 


-51.4 ± 


0.6 


-143.3 ± 


1.3 




AA O 


AGAGCTCT 


c 


-49.5 ± 


1.8 


-134.5 ± 


5.7 


-7.76 ± 0.07 


50.6 


T CTCGAGA 


d 


-50.1 ± 


0.8 


-136.7 ± 


1.8 


-7.76 ± 0.22 


50.5 


TGAGCTCA 


c 


-50.7 ± 


0.7 


-138.4 ± 


2.2 


-7.73 ± 0.04 


50.1 


ACTCGAGT 


d 


-50.3 ± 


0.7 


-137.3 ± 


1.6 


-7.72 ±0.18 


50.1 



a The top strand of each duplex is represented in the 5' to 3 1 orientation and the bottom strand 
is shown in the 3* to 5 T direction. Terminal mismatch nearest neighbors are represented in bold. 
Mismatches are underlined. T M calculated using 10* 4 M total strand concentration. 

c Thermodynamic parameters from averaging the fits of melting curves. Reported errors are 
standard deviations in the precision of the data. d Thermodynamic parameters from Tm' 1 vs. ln(C T ) 

plots. Reported errors are standard deviations in the precision propagated from the slope and 
intercept of the 1/T M vs. In C T plot. c Data from reference (19). 
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Table 1 : Thermod> namic Parameters for Hairpin Oligomer Association and Oligomer Dupiex 
Formation. 

AH 0 AS 0 AG 0 ;- T M 

(kcal/mol) (cal/molK.) (kcal mol) (°Ci 

Systems with Flementarv Interfaces" 

-AAGCCTT GA-ACAACG c -50.7 ± 4.1 -141.7 ± 11.3 -6.74 =: 0,:7 40 } 
^CGCGGAACT / TGTTGC (i) 

r A AGC C T T G T - TC AAC G c -52.8 ± 4.2 -146.3 ± 11.7 -7.43 = 0.30 53 1 
^•CGCGGAACA/AGTTGC 



GCAACT - TGTTCCGA A-v c -63.5 ± 5.1 -179.2 ± 14.3 -7.99 ± 0.32 53.2 
CGTTGA/ACAAGGCCC^ 



AAGCCTT GA - TCAACG c -53.6 ± 4.3 -149.3 ± 11.9 -7.34 ± 0.29 52.3 
CGCGGAACT/AGTTGC <m) 



-AAGC C TT GT - ACAACG c -45.1 ± 3.6 -124.8 ± 10.0 -6.42 ± 0.26 48.4 
^•CGCGGAAC A/TGTTGC 

A AGC C T T G C - AC AAC G c -46.1 ± 3.7 -128.5 ± 10.3 -6.26 ± 0.25 46.9 
CGCGGAA C G / TGTTGC 

A AGC C T T G T - GCAAC G e -52.2 ± 4.2 -144.4 ± 11.5 -7.39 ± 0.30 53.1 
CGCGGAACA/CGTTGC 

AAGCCTT GG- TCAACG e -53.6 ± 4.3 -148.2 ± 11.9 -7.67 ± 0.31 54.4 
CGCGGAA C C / AGTTGC 

AAGCCTT GA-CCAACG e -51.3 ± 4.1 -140.2 ± 11.2 -7.81 ± 0.31 56.2 
CGCGGAACT / GGTTGC 



c 



c 



c 



c 
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Table 1: Continued. 



AH 0 AS 0 AG° 3 , t m 

(kcal/mol) (cal / mol K) (kcal mol) (°C> 

-A AGC C TT GC - TCAACG c -46.1 * 3.7 -126.3 * 10.1 -6.90 * 0.:8 51." 
^CGCGGAACG/AGTTGC 



.AAGCCTTGA-GCAACG c -48.2 ± 3.9 -131.7 ± 10.5 -7.32 ± 0.29 53.9 
^- C GCGGAA C T / CGTTGC 



GC A AC A - G G TTCCGAA-v c -51.9 * 4.2 -147.0 * 11.8 -6.34 * 0.25 46.3 
CG TTGT / C CAAGGCCC^ 

-AAGCCTTGG-ACAACG 6 -50.4 * 4.0 -139.5 * 11.2 -7.10 * 0.28 51.7 
^- C GCGGAA CC/TGTTGC 

AAGCCTT GT - CCAACG 6 -54.2 ± 4.3 -147.9 * 11.8 -8.29 * 0.33 58.3 
C GCGGAA C A / GGTTGC 

AAGCCTT GC-GCAACG 6 -47.6 * 3.8 -130.9 * 10.5 -6.96 * 0.28 51.6 
CGCGGAACG/CGTTGC 

A AGC C T T GG - CC AAC G e -52.1 * 4.2 -140.1 * 11.2 -8.67 * 0.35 61.9 
C GCGGAA CC /GGTTGC 

AAGCCTTGC-CCAACG '-53.3 * 4.3 -146.5 * 11.7 -7.92 * 0.32 56.1 
C GCGGAA CG /GGTTGC 

AAGCCTT GG-GCAACG e -49.3 * 3.9 -135.2 * 10.8 -7.35 * 0.29 53.8 



c 



c 



c 



c 



c 



CGCGGAACC/CGTTGC 
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Table 1: Continued. 



(kcal / mol) 

Svstwtk; with Dangling Ends at the Interface" 

AAGCCTT GC - GCAACG c -44.4 ± 3.6 



AS 0 



AG°3-r t m 



c 



c 



c 



c 



c 



c 



c 



c 



c 



(cal / mol K.) (kcal ; mol) 



-121.4 * 9.7 -6.76 * 0.27 51 1 



CGCGGAA CG / CGTTGC 
A 



AAGCCTT GC-GCAACG e -48.0 * 3.8 -131.8 ± 10.5 -7.16 * 0.29 52.9 



CGCGGAACG/CGTTGC 

A 



AAGCCTTGC-GCAACG c -46.3 ± 3.7 -127.4 * 10.2 -6.83 ± 0.27 51.0 



CGCGGAACG/CGTTGC 
T 



AAGCCTTGC-GCAACG e -49.0 * 3.9 -136.6 * 10.9 -6.59 ± 0.26 48.7 



CGCGGAACG/CGTTGC 

T 



AAGCCTTGC-GCAACG c -37.6 ± 3.0 -102.0*8.2 -5.91* 0.24 46.4 



CGCGGAACG/CGTTGC 
A A 

AAGCCTTGC-GCAACG c -36.2 ± 2.9 -97.3 ± 7.8 -6.03 ± 0.24 47.8 
CGCGGAACG/CGTTGC 
T T 

AAGCCTTGC-GCAACG e -44.0 ± 3.5 -123.5 ± 9.9 -5.67 ± 0.23 43.0 
CGCGGAACG/CGTTGC 
A T 

AAGCCTTGC-GCAACG e -43.6 ± 3.5 -119.5 ± 9.6 -6.53 ± 0.26 49.6 
CGCGGAACG/CGTTGC 
T A 

AAGCCTTGG-TCAACG c -47.2 ± 3.8 -129.8 * 10.4 -6.90 * 0.28 51.3 
CGCGGAACC/AGTTGC 

A 
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Table 1: Continued. 



AH 0 AS 0 AG° 37 T M 

(kcal/mol) (cal/molK) (kcal / mol) (°C) 

Systems with Extra Central Nucleotide at the Interface" 

r A A GC C T T G C AGC AAC G c -44.4 ± 3.6 -122.2 ± 9.8 -6.50 * 0.26 49.2 
^ C GCGGAA C G / CGTTGC 

AAGCCTTGTACCAACG s -45.0 * 3.6 -124.6 ± 10.0 -6.32 ± 0.25 47.7 



Qgcggaaca/ggttgc 

Oligomers 



TCAACG e -38.5 ± 3.1 -108.3 ± 8.7 -4.94 ± 0.20 38.0 
AGTTGC 

AC AAC G c -36.1 ±2.9 -99.3 ± 7.9 -5.32 ± 0.21 41.4 
TGTTGC 

GCAACG c -42.5 ± 3.4 -117.1 ± 9.4 -6.21 ± 0.25 47.5 
CGTTGC 

CC AAC G e -38.6 ± 3.1 -106.2 ± 8.5 -5.69 ± 0.23 44.2 
GGTTGC 

TGTTGC c -37.1 ± 3.0 -101.1 ± 8.1 -5.79 ± 0.23 45.3 
AC AAC G 

AGTTGC c -37.0 ± 3.0 -101.0 ± 8.1 -5.70 ± 0.23 44.2 
TCAACG 
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4 T M calculated using 10"* total strand concentration. 

b The top strand of each system is conventionally represented in the 5* to 3' orientation. Nucleotides 
involved in coaxial stacking interfaces are represented in bold. 

c Parameters obtained by averaging the results of melt fit and TM" vs. ln(C T /4) plot methods. 
Errors are estimated to be 8% for AH° and AS 0 and 4% for AG° 37 . 
(i). do, (no M e j t j n g curves for these systems are shown in Figure 2. 



64 



Table 2: Thermodynamic Parameters for Coaxial Stacking*. 



AH°(coaxial stacking) 
(kcal / mol) 

Flffflrmar' interfaces* 

GA-AC -14-6 ± 5.0 



AS°(coaxial stacking) AG° 37 (coaxial sucking) 
(cal / mol K) fccal / mol) 



CT / TG 

GT - TC 

C A / AG 

CT - TG 
GA / AC 

GA - TC 
CT / AG 

GT - AC 
CA / TG 

GC - AC 

CG / TG 

GT - GC 

CA / CG 

GG - T C 

CC / AG 

GA - CC 
CT / GG 

GC - TC 
CG / AG 

GA - GC 

CT / CG 

CA - GG 
GT / CC 

GG - AC 
CC / TG 



-14.3 ± 5.2 
-26.6 ± 5.9 
-15.1 ± 5.3 
.9.0 ± 4.6 
-10.0 ± 4.7 
-9.6 * 5.4 
-15.1 * 5.3 
-12.7 ± 5.1 
-7.6 * 4.8 
-5.6 ± 5.1 
-14.8 ± 5.1 
.14.2 * 5.0 



-42.4 ± 13.8 
-38.0 ± 14.6 
-78.2 ± 16.5 
-41.0 ± 14.8 
-25.5 * 12.8 
-29.2 ± 13.0 
-27.3 ± 14.9 
-39.9 ± 14.7 
-34.0 ± 14.1 
-18.0 ± 13.3 
-14.6 ± 14.1 
-45.9 * 14.3 
-40.2 ± 13.7 



-1.42 * 0.34 
-2.49 ± 0.36 
-2.29 ± 0.39 
-2.40 ± 0.35 
-1.10 ± 0.33 
-0.94 ± 0.33 
-1.18 ± 0.39 
-2.73 * 0.36 
-2.12 ± 0.39 
-1.97 ± 0.34 
-1.11 ± 0.38 
-0.56 ± 0.34 
-1.78 ± 0.35 
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Table 2: Continued. 



AH°(coaxial stacking) AS°(coaxial stacking) AG° 37 (coaxial stacking) 
(kcal / mol) (cal / mol K.) (kcal / mol) 

GT-CC -15.6 ± 5.3 -41.8 ± 14.6 -2.61 ± 0.40 
CA / GG 

GC -GC -5.0 ± 5.1 -13.8 ± 14.1 -0.75 ± 0.37 
CG / CG 

GG-CC -13.5 ± 5.2 -33.9 ± 14.1 -2.98 ± 0.41 
CC / GG 

GC-CC -14.7 ± 5.3 -40.3 ± 14.5 -2.23 ± 0.39 
CG / GG 

GG-GC -6.8 ± 5.2 -18.1 ± 14.3 -1.14 ± 0.38 
CC / CG 

Interfaces with Dangling Ends b 

GC-GC -1.9 ± 4.9 -4.3 ± 13.5 -0.55 ± 0.37 

CG / CG 
A 

GC-GC -5.5 ± 5.1 -14.7 ± 14.1 -0.95 ± 0.38 

CG / CG 
A 

GC-GC -3.8 ± 5.0 -10.3 ± 13.8 -0.62 ± 0.37 

CG / CG 
T 

GC-GC -6.4 ± 5.2 -19.5 ± 14.4 -0.38 ± 0.36 

CG / CG 
T 

GC-GC 5.0 ± 4.5 15.1 ± 12.4 0.30 ± 0.34 

CG / CG 
A A 
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Table 2: Continued. 



AH°(coaxial stacking) AS°(coaxial stacking) AG° 37 (coaxial stacking) 
(kcal / mol) (cal / mol K) (kcai / mol) 

GC-GC 6.3 ± 4.5 19.8 ± 12.2 0.18 ± 0.35 

CG / CG 
T T 

GC-GC -1.4 ± 4.9 -6.4 ± 13.6 0.55 ± 0.34 

CG / CG 
A T 

GC-GC -1.1 ± 4.9 -2.4 ± 13.4 -0.32 ± 0.36 

CG / CG 
T A 

A 

GG-TC -8.6 ± 4.9 -21.5 ± 13.5 -1.96 ± 0.34 

CC / AG 
A 

Interface with Fvtra Central Nucleotide'* 

GCAGC -1.9 ± 4.9 -5.1 ± 13.5 -0.29 ± 0.36 
CG / CG 

GTACC -6.4 ± 4.7 -18.5 ± 13.1 -0.64 ± 0.34 
CA / GG 

* These parameters and their corresponding errors are deduced from Table 1 as described in the text. 
b The top strand of each duplex is conventionally represented in the 5' to 3' orientation. Nucleotides 
involved in coaxial stacking interfaces are represented in bold. 
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Table SI: Extinction coefficients of hairpins at 25 °C 



C 



A AGC c ttggtcaacg 

C GCGGAACC 



gyperimental' 
(L mol' 1 cm' 1 ) 
188847 



calculated 
(L mol' 1 cm" 1 ) 
192310 



c 



AAGC cttgatcaacg 

C GCGGAACT 



188718 



195950 



c 



AAGCCTTGCTCAACG 
C GCGGAACG 



186139 



191610 



c 



AAGC cttgttcaacg 

CGCGGAACA 



191071 



195810 



c 



AAGCCTTGAACAACG 
C GCGGAACT 



194330 



200950 



c 



AAGCCTTGTACAACG 
C GCGG AACA 



193953 



200810 



c 



AAGC C TTGGACAACG 
C GCGGAACC 



188889 



197310 



c 



AAGCCTTGCACAACG 
C GCGG AACG 



194623 



196610 



C 



A AGC C TTGAGCAACG 
C GCGG AACT 



192953 



197350 



C 



AAGC CTTGTGC AACG 
CGCGGAACA 



195968 



197630 



C 



AAGCCTTGGGCAACG 
C GCGGAACC 



195177 



193710 
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Table SI: Continued 



c 



A AGCCTTGCGCAACG 
CGCGGAACG 



190944 



193010 



C 



AAGC c ttgaccaacg 

C GCGGAACT 



192663 



194350 



C 



AAGCCTTGTCCAACG 
C GCGG AACA 



193532 



194210 



C 



AAGC CTTGGCCAACG 
C GCGG A ACC 



195094 



192390 



C 



AAGCCTTGCCCAACG 
CGCGGAACG 



192864 



190010 



GC AACA - G TTCC AA-s 
CCAAGGCCC^ 



200806 



196010 



GC AACT - T TTCC AA>, 
AC AAGGCC C^ 



206650 



192810 



C 



AAGCCTTGCAGCAACG 
CGCGGAACG 



202688 



206510 



C 



AAGCCTTGTACCAACG 
C GCGG AACA 



204450 



208010 



C 



AAGC CTTGCGCAACG 
C GCGGAACG 
A 



191282 



193010 



C 



AAGCCTTGCGCAACG 
C GCGGAACG 
T 



198343 



193010 
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4 Calculated with Equation 3. 
b Calculated with Equation 4. 



O 

TO 
S3 

m 
m 

4 
« : 

O 
01 
O 

o 
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Table S2. Thermodynamic Parameters for Hairpin Oligomer Association and Oligomer Dupiex 
Formation. 



AH 0 



V 



c 
c 
c 



Eiementarv interfaces' 

AAGCCTTGA- ACAACG 0> 
CGCGGAACT/TGTTGC 

AAGCCTTGT-TCAACG(«) 
CGCGGAACA / AGTTGC 

AAGCCTTGA- TCAACG m 
CGCGGAACT /AGTTGC 



GCAACT-TGTTCCGAAn 
CGTTGA / ACAAGGCCC-' 



C 

c 



AAGCCTTGT-ACAACG 
CGCGGAACA /TGTTGC 

AAGCCTTGC-ACAACG 
CGCGGAACG/TGTTGC 



-AAGCCTTGT -GCAACG 
^CGCGGAACA/CGTTGC 

-AAGCCTTGG-TCAACG 
GCGGAACC /AGTTGC 

r AAGCCTTGA-CCAACG 
^CGCGGAACT /GGTTGC 



C 

c 



A AG CCTTGC-TCAACG 
C GCGGAACG /AGTTGC 

AAGCCTTGA-GCAACG 
CGCGGAACT /CGTTGC 



(kcal / mol) 


(cal / mol K) 


(kcal / mol") 


( C C) 


-48.1 ± 


2.5 


-133.2 


± 


8.1 


-6.8 




OA 


50. 1 


-53.3 ± 


2.7 


-150.1 


± 


8.7 


-6.7 


± 


0.0 


48 4 


-55.7 x 


4.3 


-155.8 


± 


14.0 


-7.4 


± 


0.2 


52.I 


-49.9 ± 


1.7 


-136.9 


± 


5.6 


-7.4 


± 


0.0 


54.2 


CIO _i_ 


4.0 


-143.6 


± 


12.7 


-7.3 


± 


0.1 


52.8 


-55.4 x 


2.1 


-155.0 


± 


4.9 


-7.3 


± 


0.5 


51.7 


-61. 6 ± 


2.6 


-173.1 


± 


7.9 


-8.0 


± 


0.1 


53.7 


-65.5 ± 


3.1 


-185.2 


± 


9.8 


-8.0 


± 


0.0 


52.8 


-48.0 ± 


2.3 


-134.4 


± 


7.2 


-6.4 


± 


0.1 


47.3 


-42.2 ± 


1.9 


-H5.2 


± 


6.3 


-6.5 


± 


0.1 


49.5 


-45.9 ± 


1.1 


-127.9 


± 


3.3 


-6.3 


± 


0.1 


47. 1 


-46.3 ± 


4.0 


-129.2 


± 


I3.'l 


-6.3 


± 


O.l 


46.8 


-54.7 ± 


2.7 


-152.4 


± 


9.3 


-7.4 




0.2 


52.2 


-49.7 ± 


3.4 


-136.3 


± 


ll.O 


-7.4 




0.1 


53.9 


-53.5 ± 


4.3 


-147.7 


± 


13.6 


-7.7 


± 


0.2 


54.4 


•53.8 ± 


3.4 


-148.7 


± 


111 


-7.7 


± 


0.1 


54.3 


-52.9 =fc 


2.2 


-145.4 


± 


7.3 


-7.8 


± 


0.1 


55.6 


-49.6 ± 


2.4 


-134.9 


± 


5.5 


-7.8 


± 


0.7 


56.8 


-48.6 ± 


2.8 


•134.6 


± 


9.3 


-6.9 


± 


0.2 


50.6 


-43.6 ± 


22 


-H8.1 


± 


7.3 


-6.9 


± 


0.1 


52.8 


-49.7 ± 


3.5 


-136.6 


± 


ll.4 


-7.3 


± 


0.1 


53.4 


-46.6 ± 


1.8 


-128.8 




5.7 


-7.3 


± 


0.0 


49.9 
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Table S2: Continued. 



GCAACA -GGTTCCGA A-v 
CGTTGT / CCAAGGCCC^ 



c 



c 



c 



c 



c 



AAGCCTTGG-ACAACG 
CGCGGAACC / TGTTGC 

AAGCCTTGT-CCAACG 
C GCGGAACA / GGTTGC 

AAGCCTTGC-GCAACG 
CGCGGAACG/ CGTTGC 

AAGCCTTGG-CCAACG 
CGCGGAACC/ GGTTGC 

AAGCCTTGC-CCAACG 
CGCGGAACG /GGTTGC 



C 



AAGCCTTGG-GCAACG 
CGCGGAACC /CGTTGC 

Interfaces with dangling ends* 



AH° 
(kcal / mol) 



-53.1 
-50.8 



-48.4 ± 
-52.3 ± 

-58.0 ± 
-50.3 ± 



-57.3 ± 
-46.9 ± 

-55.3 ± 
-51.4 ± 



3.8 
3.3 

3.1 
2.0 

2.1 
2.1 



-47.0 ± 3.9 
-48.2 ± 0.9 



5.4 
2.7 

2.6 
1.5 



-45.9 ± 1.0 
-52.6 ± 3.5 



AS 0 
(cal / mol K) 



-150.8 
■143.3 

■133.1 
■145.9 

■160.2 
■135.7 



■129.0 ± 
•132.8 ± 



• 156.4 
-123.9 

-152.7 
•140.2 

■124.4 
-146.0 



13.0 
10.9 

10.4 
6.5 

6.4 
4.7 

12.9 
2.9 

17.0 
5.8 

8.1 
3.4 

3.5 
11.3 



AG° 3 - 
(kcal mol) , 



T 



-6.3 - 0.2 

-6.4 = 0.1 

-7.1 ± 0.1 

-7.1 ± 0.0 



-8.4 
-8.2 

-7.0 
-7.0 



-7.9 
-7.9 



-8.8 ± 
-8.5 ± 



-7.3 ± 
-7.4 ± 



0.2 
0.6 

0.1 
0.0 

0.3 
0.9 

0.1 
0.5 

0.1 
0.1 



4 

4< 



3 
5 

5 
5 

6 
6 

5 
5 

5 
5 



c 



AAGCCTTGC-GCAACG « -45.2 ± 2.9 -123.8 ± 9.4 
CGCGGAACG/CGTTGC d -43.7 ± 3.4 -119.1*11.1 
A 



-6.8 ± 0.2 5 
-6.8 ± 0.1 51 



C 



AAGCCTTGC-GCAACG « -46.7 ± 3.7 



CGCGGAACG/CGTTGC 
T 



•46.0 ± 3.3 



-128.7 ±11.8 -6.8 ± 0.1 50 
•126.2 ± 7.5 -6.8 ±1.0 5 



-AAGCCTTGC - GCAACG 
^CGCGGAACG/CGTTGC 

A 



< -48.6 ± 



3.0 

-47.4 ± 3.2 



•133.7 ± 9.4 
•129.8 ± 10.5 



-7.2 ± 0.1 
-7.2 ± 0.1 



72 




Table S2: Continued. 







AH° 




AS 0 


AG 0 


3 " 


T 1 

* M 






(kcal / mol) 


(cal / mol K) 


(kcal / 


mol) 


(°C) 


,-AAGCCTTGC-GCAACG 
^CGCGGAACG / CGTTGC 

T 


C 

d 


-46.3 
-51.6 


± 
± 


2.6 
3.7 


-127.9 
-145.4 


± 


8.7 
12.1 


-6.6 
-6.5 




0.2 
0.1 


43 6 
42.3 


-AAGCCTTGG- TCAACG 
^-C GCGGAACC / AGTTGC 

A 


C 

d 


-49.2 
-45.1 


± 
± 


3.5 
2.9 


-136.6 
-123.1 


± 
± 


11.0 
9.4 


-6.9 
-6.9 


± 


0.2 
0.1 


44.8 
45.7 


--AAGCCTTGC - GCAACG 
^CGCGGAACG/ CGTTGC 

▲ A 

A A 


c 
d 


-39.7 
-35.5 


± 
± 


4.6 
2.9 


-109.1 
-95.0 


± 
± 


15.2 
5.9 


-5.8 
-6.0 


± 
± 


0.2 
1.1 


38.2 
39.8 


-AAGCCTTGC-GCAACG 
^•CGCGGAACG / CGTTGC 

A T 


c 
d 


-44.7 
-43.2 


± 
± 


3.3 
3.7 


-126.0 
-121.0 


± 
± 


11.0 
12.2 


-5.6 
-5.7 


± 
± 


0.2 
0.2 


36.6 
37.2 


r AAGCCTTGC-GCAACO 
^- C GCGGAACG / CGTTGC 

1 1 


c 
d 


-38.5 
-33.9 


± 
± 


3.0 
1.3 


-105.1 
-89.6 


± 
± 


10.0 
4.5 


-6.0 
-6.1 


± 
± 


0.2 
0.1 


39.3 
41.0 


AAGCCTTGCftfAArrf 
^C GCGGAACG /CGTTGC 
T A 


c 
d 


-43.3 
-43.9 


± 


9.0 
1.0 


-118.5 
-120.5 


± 
± 


30.1 
3.4 


-6.5 
-6.5 


± 
± 


0.0 


43.1 
43.2 


Interfaces with extra central nucleotide 
























-AAGCCTTGCAGCAACQ 
^C GCGGAACG /CGTTGC 


e 
d 


-43.8 
-45.0 


± 
± 


9.4 
3.5 


-1202 
•124.1 


± 
± 


31.0 
11.5 


-6.5 
-6.5 


± 
± 


A ^ 
V.Z 

0.1 


42.7 


-AAGCCTTGTACCAACG 
^CGCGGAACA/GGTTGC 


c 
d 


-45.0 ± 
-45.0 ± 


2.8 
4.0 


-124.5 
-124.8 


± 
± 


9.1 
13.4 


-6.3 
-6.3 


± 


0.2 
0.2 


41.6 
41.4 
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Table S2: Continued. 



AH 0 AS 0 AG° 37 V 

(kcal / mol) (cal / mol K) (kcal / mol) { *C) 

Oligomers 



T C A A C G 


<= -40.6 




2.5 


-115.4 


± 


9.1 


-4.8 


X 


0.3 


30.5 


AG T T G C 


d -36.5 


± 


1.9 


-101.2 


± 


6.3 


-5.1 


± 


0.1 


31.9 


A C A A C G- 


« -38.3 


± 


2.0 


-106.8 


± 


6.5 


-5.2 




0.2 


jj.4 


T G T T G C 


-33.9 




1.6 


Al ft 

-91.8 


± 


5.3 


-5.4 




0.1 


34.6 ' 


GC A A C G 


« -43.9 


± 


2.4 


-121.6 


± 


7.4 


-6.2 


± 


0.1 


40.7 


CGT T G C 


4 -41.1 


± 


2.1 


-112.6 




6.9 


-6.2 


± 


0.1 


41.2 


CCAACG 


c -41.0 


± 


2.6 


-114.1 


± 


8.9 


-5.6 


± 


0.2 


36.4 


GGT T G C 


4 -36.2 


± 


0.9 


-98.2 


± 


1.8 


-5.8 


± 


0.3 


37.8 


TGT T GC 


« -38-8 


± 


3.8 


-106.6 


± 


12.1 


-5.7 


± 


0.1 


37.5 


AC A A C G 


4 -35.5 


± 


2.8 


-95.7 


± 


9.2 


-5.8 


± 


0.1 


38.4 


AGT T G C 


« -37.1 


± 


2.9 


-101.2 


± 


9.4 


-5.7 


± 


0.1 


36.9 


TC A AC G 


d -36.9 


± 


2.4 


-100.8 




8.0 


-5.7 


± 


0.1 


36.9 



1 T M calculated for 4x10"* total strand concentration. 

b The top strand of each system is conventionally represented in the 5' to 3* orientation. Nucleotides 
involved in coaxial stacking interfaces are represented in bold. 

c Parameters obtained from averaging fits of melting curves. Reported errors are standard deviations 
in the precision of the fitted data. 

d Parameters obtained from T M vs. ln(Cj/4) plots. Reported errors are standard deviations in the 
precision propagated from the slope and intercept of the 1/T M vs. In (Cj/4) plot. 
(0. <»). (iii) irT ^ vs i^Ct/4) plots for these systems are shown in Figure SI . 
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